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A Technique for Real Time Simulation 
of a Rigid Body Problem 


1. Introduction. With the advent of large-scale high-speed digital com- 
puters, there arises the question of their possible use in the solution of 
problems in “real time,” i.e., in conjunction with instruments receiving and 
responding to stimuli from the external environment. The criteria for satis- 
factory operation in such real-time service are different from those generally 
encountered. First, the computer must calculate fast enough to respond to 
the fastest changing stimuli; speed in solving the associated equations is in 
fact the primary requirement. Again, because it is impossible accurately to 
predict stimuli variations, either the response must be delayed or a “‘reason- 
able’ amount of inaccuracy must be tolerated or both. Third, because of the 
limited range of number representation in the digital computer (or, for that 
matter, in any computer), all parameters must remain bounded no matter 
when or where they occur. Finally, the computations must be stable. 

This paper is devoted to a discussion of the manipulation of the dynam- 
ical equations of a rotating rigid body with fixed center of gravity.' The 
techniques here adapted will be-seen to sacrifice accuracy to speed of 
computation, but in such a way as to hold inaccuracies reasonably low; at 
the same time, a singularity generally encountered will be circumvented. 
Stability will not be discussed. 

2. Manipulation of Eulerian angles. The equations of motion of a 
rotating rigid body are usually simplified by using a set of axes attached to 
the rigid body. The simplification results from the fact that this choice of 
axes causes the moments and products of inertia to become constants. In 
addition, the products of inertia can be made to vanish by properly orienting 
these axes with respect to the rigid body. 

The orientation of the rotating axes OX, OY, and OZ relative to a set of 
fixed axes OX», OYo, and OZ» may be defined by the Eulerian angles y, 0, ¢. 
For convenience assume OZ) downward, OX» North, and OY, East. wy is 
the angle measured about OZ, from the XZ» plane to the (vertical) plane 
containing OZ» and OX;; it is positive when OX is toward the East. @ is the 
angle between OX and the XoY» plane; it is positive when OX is above the 
XoYo plane. It is measured about an axis in the Xo Yo plane. Call this axis 
OY’. dis the angle measured about the OX axis from OY’ to OY; it is positive 
for clockwise rotation as -riewed from the origin along the positive OX axis. 

If the angular velocity of the moving axes has components , g, and r 
about the moving axes OX, OY, and OZ, respectively, then the following 
three differential equations show the relations among the angles and the 
components of angular velocity. 


6 = qcos — rsin ¢, 
1 
= (P 0088 + sin + 7 sin 8 cos ¢), 
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The presence of cos @ in the denominator makes two of these equations 
indeterminate when @ is +90°. This makes the equations difficult to use 
especially where point by point numerical computations are being performed. 
One common method is to limit the variable @ say to the range +85°. This is 
avoiding the problem rather than solving it. 

This indeterminism is a product of the mathematical approach to a 
physical problem. A different choice of axes for example would present the 
same difficulty at a different orientation of the rotating body. Therefore the 
way out of the difficulty required a new mathematical point of view. 

3. Derivation of direction cosine equations. When the equations of 
motion contain the Eulerian angles they appear in combinations equal to 
one or more of the direction cosines. That is, the Eulerian angles can be 
eliminated from the mechanical equations of motion by introducing the 
direction cosines as parameters. The kinematic differential equations for the 
direction cosines are quite simple and have no indeterminate points. This 
is the principal gain. 

The new differential equations can be derived by the following method. 
The unit vector along the Xo axis is related to the unit vectors along the 
moving X, Y, Z axes by 


to = th jm, 
Differentiate this with respect to time 
di dj dk 
0= + il, + + jm: + amt kn. 

If the moving axes are rotating with an angular velocity ip + jg + kr 

Substituting these values for the derivatives we obtain 


i(l, — + qm) + — pm + rh) + R(t — gh + pm) = 0, 
which requires that 


i, = rm, — qn, 
my, = pn rly, 
ny = ql pm,. 


A similar treatment of j) and ko yields similar equations with subscripts 2 
and 3, respectively. 

4. Identities. Since there were three independent differential equations 
involving the Eulerian angles only three of these nine new differential 
equations can be expected to be independent of the 21 identities that obtain 
for direction cosines. These 21 identities include six expressions for normality 


(1) l? + m? + n? = 1, etc., 
six for orthogonality 
(2) Lyle + mymz + nyn2 = 0, etc. 


anc 
©) 
not 
sar 
is 
ust 
en 
tio 
| 
ab 
co 
th 
he 
rc 
ra A 
| 
t 


a2 


ons 
tial 
ain 


lity 


REAL TIME SIMULATION OF A RIGID BODY PROBLEM 75 


and nine of the form 
(3) l, = mans — etc. 


To choose three independent equations, choose equations all of which do 
not have the same subscripts nor do all involve the time derivative of the 
same letter. For example 

1 = 7m, — qn, 

i, = — qn2, 

= pnz — re 

is a set of independent equations which, with six of the identities, can be 
used :o obtain all nine direction cosines. If desired, more than three differ- 
ential equations and fewer than six identities can be used. For some applica- 
tions, such as analogue computers, the differential equations are more suit- 
able than the identities. Even for digital computers the differential equations 
may be preferred to identities of the form of equations (1) and (2). 

5. Computation Procedure. Before discussing the plan of computation, 
it will be well to discuss a fundamental change that must often be made. 
Since there is no unity in many modern digital computers, this number 
cannot be held in the machine. Indeed, if a direction cosine were computed 
correctly to be unity, it would appear in such a machine as zero. To remove 
this difficulty, pseudo direction cosines can be introduced, which are M 
times the corresponding true direction cosines. For example, with M = 1/2, 
the identities are of the form 


13? + + = (6 of these), 
Isle + mym2 + nz = 0 (6 of these), 
l, = — mgnz) (9 of these). 

The first step in the proposed computation procedure is the use of three 
differential equations to compute /3;, m3, m3. Since both curtailment and 
round-off errors will gradually cause these parameters to fail to obey the 
normality equations, a second step is required to correct this failure. 

Let the computed values be /;’, m;’, m3’. Define « by the equation: 

1s’? + ms" + n;" = M*(1 + ©). 
A normalized si of direction cosines is therefore 
2 
nz = i+ ee 
By using two terms in the binomial expansion of (1 + «)-+, 


1 
= (1 = (5 = ‘), similarly ms, 


thus avoiding the square root computation. 
1,', ms’, m2’ may now be computed from the appropriate three differential 
equations. Consider the two vectors, 
r3 = ils + jms + 
re’ = il,’ + jm! + 
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The vector r3 has magnitude M. The vector 7,’ should, but in general 
will not, have magnitude M and be perpendicular to 73. Let re’ be turned 
through as small an angle as possible to make it perpendicular to r3. To do 
this it is necessary only to look for r2 in the plane of re’ and 73, i.e., 


= (i a)r2’ + 


The scalars a and 6 must at the same time make the magnitude of 7. 
equal to M. 
Define ¢; and e2 by the equations 


+- me’? + ne” = M*(1 + 


If the errors €; and €2 are small, the corrections a and 6 will be small; then 
it will be sufficiently exact to neglect second degree terms in a, b, «1, €:. 
Requiring r2 to have magnitude M leads to 


while the fact that rz is normal to r; yields 
b = — 


Thus, neglecting second order terms, the corrected values may now be 
written 


l, = 


| 
| 
| 


m2 = 


= (1 — N3€2. 


The correction involving €; adjusts the magnitude while the correction in- 
volving €2 corrects the direction of r2. The results obtained show that to the 
first order either correction may be applied first or they may both be 
applied together. 

There is a second method of orthogonalizing, which is preferable because 
it is quicker. In this method the computed values of two of 12, m2, m2 are 
assumed correct and the third is adjusted to make r2 perpendicular to 73. 
If m3 is greater than either /; or ms, for example, only /; and mz are computed 
by integrating their differential equations; 2 is then determined from the 
orthogonality condition 
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Note that if m2’ had been computed from the differential equation, then 


lols + moms + = 
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This amounts to correcting m2’ by the term 


€3;M? 


Obviously if J,’ had been corrected instead of m2’ the correction would be 


— 
ls 


To make the correction small, the largest of /;, ms, m3 is made to appear in 
the denominator. This also prevents division by zero, which might occur if 
m2 were corrected under all circumstances. 

After the direction of rz has been established perpendicular to rz, its 
magnitude must be adjusted. The procedure is to alter all components pro- 
portionately as for r; above. Clearly this does not destroy orthogonality. 
l,, my, m, can now be computed from equation (3) 


mn; — msn 


6. Conclusions. The equations of motion of a rotating rigid body can be 
formulated directly in terms of the direction cosines. The direction cosines 
can be determined as described above in a manner which is thought to be 
ideally suited to digital computation. No trigonometric functions of the 
Eulerian angles appear, obviating the necessity of referring to a function 
table or to a series expansion to obtain the solutions of the equations of 
motion. li addition, all quantities are finite for all orientations of the 
rigid body. 

This work was done in connection with a study performed by the Uni- 
versity of Pennsylvania, Moore School of Electrical Engineering, under 
Contract N6Onr 24913 sponsored by the Office of Naval Research, Special 
Devices Center, Port Washington, New York. 


H. J. Gray, Jr. 
M. RUBINOFF 
H. SoHON 
University of Pennsylvania 
Moore School of Electrical Engineering 


1 LetGu PaGE, Introduction to Theoretical Physics, Second Edition—Fourteenth Printing, 
D. Van Nostrand Company, Inc., New York, N. Y., Introduction, Chap. I, II; 1951. 


A Solution of Simultaneous Linear Equations 
and Matrix Inversion with High Speed 
Computing Devices 


In solving several systems of simultaneous linear equations, in which 
the coefficients of the unknowns are the same in all the systems but in 
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which the constant terms vary from system to system, it is advantageous 
to invert the matrix formed by the coefficients of the unknowns. This in- 
verted matrix, common to all the systems, is then multiplied by the vector 
formed by each set of constant terms, and a set of solutions is thus obtained 
for each of the above systems,! or simultaneous solutions of the systems ob- 
tained at the time of inversion. The inversion of a matrix, especially if it is 
of a high order, is a laborious and time-consuming process. From time to 
time new methods of inversion have been developed and modified to suit 
the computation equipment available until it seems that a saturation point 
in simplicity of application, accuracy of results, and speed of computation 
has been reached, and further improvements will have to come from the 
development of the computing machines themselves or from their applica- 
tion. For information on these topics the reader is referred to items 2-12 of 
the bibliography. 

The following is a modification of an algebraic treatment of the matrix 
by the direct method which will reduce the labor and time necessary for 
inversion. This method has time and labor saving characteristics applicable 
to computations with low and high memory computing devices, making 
the inversion of high order matrices less forbidding. 

In addition to the customary definitions of inversion, pivotal row, lead- 
ing element, etc., which are here taken for granted, we will define prepivotal 
and transpivotal rows as the rows respectively preceding and following the 
pivotal row. Consider also the original matrix as extended to the right by 
means of the identity matrix and a check column (other augmentations of 
the original matrix may be included to suit the needed results® and a check 
row may be added for ease in locating the source of spurious results). 

In a simplified method" of inversion the elements of this extended matrix 
are operated upon by certain reduction operations and an intermediate 
matrix is obtained, the elements of which, on a second application of the 
reduction formulas, yield another intermediate matrix, and so on until, 
after the mth reduction, the inverted matrix is obtained. The reduction 
operations usually consist of a quotient operation applied to the elements of 
the pivotal row and a binomial operation applied to the elements of the 
other rows. This reduction cycle has a definite advantage in its simplicity, 
but it involves computations with the elements of the transpivotal rows in 
all the previous reductions, in preparation of values for the time when these 
transpivotal rows become pivotal or prepivotal. Also, with automatic com- 
puting machines, all the elements of the identity matrix are involved in 
the calculations even at the time when some of them do not contribute to 
the final result. 

The method proposed in this article computes partial reductions (on 
some of the rows only) on pivotal, prepivotal, and the first transpivotal rows, 
until the memory capacity of the computing equipment is reached, at which 
point a complete reduction (on all the rows of the matrix) is computed, 
repeating this cycle, with the last obtained matrix, until final reduction 
(inversion) is reached; all this while, at a given time only the relevant con- 
tributing part of the identity matrix is used. It is to be noted that the saving 
in the identity matrix is independent of the transpivotal row economy and 
may be used with other inversion methods as well. 
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Three types of computing operations are used in each matrix reduction. 


1. Divisional or quotient operation, used on the elements of the pivotal 
row. 

2. Binomial operation, used on the elements of the prepivotal rows. 

3. Polynomial operation, used on the elements of the first transpivotal 
row for partial matrix reductions, and on all the transpivotal rows 
for complete matrix reductions. 


Let the elements of the original extended matrix be x;; (where ¢ and j 
respectively represent the row and column of the particular element, i may 
assume values from 1 to m, and j values from 1 to 2m + 1). Let the elements 
of the matrix developed after the first, second . . . and mth reductions be 
a;;, bi; ++ - mij, respectively. As usual the ith row is pivotal for the ith reduc- 
tion. The formulas for the first five reductions are: 


(1) 


= %1j/%11 


(j = 2(1)n + 1) 


by = dy — 
(2) ba; = (j = 3(1)" + 2) 
bay = X33 — — 
= Diz — 
Coz = ba; — desea; 


(3) (j = 4(1)n + 3) 
Cag = — — Xa2Cog — 
di; = Cy — 
do; = C23 — 
(4) d3; = C33 — G 5(1)n + 4) 
da; = Caj/Cas 
ds; = X53 — — — — 
= dy — dyes; 
= doz — doses; 
(j = 6(1)n + 5) 


= dy — 
= ds;/dss 


The check column should be included in the above values of j. 

Thus, on the first reduction, the pivotal row is the first row, the elements 
of which are divided and the binomial operation is applied to the elements of 
the second row. On the second reduction, the second, first, and third rows 
are pivotal, prepivotal, and the first transpivotal, respectively. In the third 
reduction, the third row is pivotal, the first and the second rows are pre- 
pivotal, and the fourth row is the first transpivotal, etc. 
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With an unlimited, or sufficiently large, memory machine the process is 
continued until the mth reduction, when the mth row is pivotal, the rows less 
than are prepivotal and there are no transpivotal rows. In applying the 
above formulas to the elements of the original (extended) and intermediate 
matrices, only columns of the identity matrix up to and including the unit 
element of the corresponding pivotal row need be included. 

With limited memory machines the process is continued until the memory 
capacity is reached, at which point the quotient operation is applied to the 
elements of the pivotal row, the binomial operation is used on the elements 
of the prepivotal rows and the polynomial operation is computed for the 
elements of all the transpivotal rows. Thus, if in the above equations the 
fourth reduction had been complete, the general equation for the transpivotal 
row elements would have been: 


dig = — — — — (¢ = 5(1)n) 


With a minimum memory machine all the reductions are computed as 
usual with a quotient operation on the elements of the pivotal row and a 
binomial operation on the elements of all the other rows. Unlimited, limited, 
and minimum memory machines are left undefined except in their use to 
classify the three broad applications of the method as described above. 
The polynomial expression for the transpivotal rows, which is easily 
obtained algebraically from previous binomial reductions by a rearrange- 
ment of terms and by a substitution of pivotal and prepivotal results of the 
same column, permits partial reductions of the matrix due to the fact that 
it depends only on the pivotal and prepivotal results of the same column 
and on the leading elements of the original extended matrix for the corre- 
sponding row. 

It is to be noticed that in the above formulas we must distinguish be- 
tween two types of leading elements: (a) derived leading elements for pivotal 
and prepivotal rows and taken from the previous matrix reduction, and (b) 
original leading elements for transpivotal rows and taken from the original 
extended matrix (or from an intermediate matrix obtained in a previous 
complete reduction). For a given transpivotal row, there are as many 
original leading elements as the number of the reduction. Varying degrees 
of continuity within the cycle itself may be obtained according to the 
storing capacity, programming facilities, available means for selection of 
properly coded rows and columns, and other characteristics of the available 
equipment. 

Between the two extremes of no memory and unlimited memory there 
is an optimum program design which is approached by a judicious choice 
of continuity of operation using the above minimum memory requirements 
to meet the available memory capacity. 

We have successfully and economically applied the above method to the 
inversion of 20 by 20 matrices using a full capacity 602-A.IBM calculating 
punch while using words of 12 digits throughout in order to take care of 
creeping decimal errors from the right and overflowing on the left. Three 
partial and one complete reductions per cycle were obtained. 

With reduced numbers, when possible, greater savings may be achieved, 
but they are not advisable when handling high order matrices. Even when 
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large integral values may be avoided by preliminary or intermediate treat- 
ment of the matrix, the creeping decimal error is ever present and may be 
met only with otherwise unnecessary decimals or with iterative methods.* 
The time consumed by this method and the time by the method of all 
complete reductions were in the ratio of 1 to 3 when using the same equip- 
ment in both cases and with the same number handicap, by the same 
operators. 
GABRIEL G. BEJARANO 
Bruce R. ROSENBLATT 
California Research Corporation 
Richmond, California 


1A. C. AITKEN, Determinants and Matrices. Edinburgh, 5th ed., 1948. 
2P.S. Dwyer, Linear Computations. New York, 1951. 
ait H. ee, “Some new methods in matrix calculation,” Annals Math. Stat., v. 14, 
1943, p. 1-34. 
4P. D. Crout, “A short method for evaluating determinants and solving systems of 
linear equations with real or complex coefficients,” AIEE, Trans., v. 60, 1941, p. 1235-1241. 
5G. W. Perris III, “Matrix inversion and solution of simultaneous linear algebraic 
equations with the IBM type 604 - Calculating Punch,” IBM Computation 
Seminar, August 1951, Proc., p. 105-11 
$j. CHANCELLOR, J. W. SHELDON, & Sc L. Tatum, “The solution of simultaneous linear 
equations using the IBM Card- Programmed Electronic Calculator,” IBM Industrial Com- 
putation Seminar, September 1950, Proc., p. 57-61. 
. BELL, “Punched card techniques for the solution of simultaneous equations 
and other matrix operations,” IBM Scientific Computation Forum, 1948, Proc., p. 28-31. 
8 J. von NEUMANN & H. H. Gotpstine, “ Numerical inverting of matrices of hig ‘order, v 
Amer. Math. Soc., Bull., v. 53, 1947, p. 1021-1099. 
* A. OPLER, “Monte Carlo matrix calculation with punched card machines,” MTAC, 
v. 5, 1951, p. 115-120. 
10 J. SHERMAN, “Computations of inverse matrices by means of IBM machines,” IBM 
Technical Newsletter, no. 3, 1951, p. 22-25. 
11 J. Lowe, “Solution of simultaneous linear algebraic equations using the IBM type 
604 Electronic Calculating Punch,” IBM Computation Seminar, Dec. 1949, Proc., p. 54- 
2 F. M. Verzun, “The solution of simultaneous linear equations with the aid of the 602 
Calculating Punch,” MTAC, v. 3, p. 453-462. 


Error Bounds on Approximate Solutions 
to Systems of Linear Algebraic Equations 


1. The Abmatrix. Consider the real matrix A = (a;;), i = 1(1)m, j = 
1(1)m. Define a(A), read “the abmatrix of A,” to be 


a(A) = 


Thus, a(A) is the matrix having each element equal to the absolute value 
of the corresponding element of A. In particular, for a scalar s, 


a(s) =|s|. 
If A = and B = (b;;), = 1(1)m, = 1(1)n, define 


a(A) < a(B) 
to mean 
< | bs; | for all i, j 


Thus, if a(A) <¢ a(B) and a(B) < a(C), then a(A) <¢ a(C). 
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From these definitions, it follows that if the operations A + B and AB 
are permissible, then 


(1) a(A + B) ¢ a(A) + a(B), 
and 
(2) a(AB) < a(A)a(B). 


Extensions to and combinations of (1) and (2) are indicated. In particular, 
if A is a square matrix and is a positive integer, 


a(A?) < [e(A)}. 


Subsequently, [a(A)]? will be written a(A). 

This simple notion of the abmatrix is readily applicable to examining 
the convergence of a class of iterative methods for solving systems of linear 
algebraic equations and to determining error bounds for approximate solu- 
tions of these linear systems; however, in this note, only the latter will be 
examined. For their application, the use of an approximate inverse is 
required. 

2. The Error in the Approximate Solution. Let the linear system under 
consideration be 


At =1, 


where A is non-singular and of order n. Let G be an approximate inverse of 
A, the approximate solution being 


x = Gn. 

Then, 

(3) (I — + (I — 

h=0 

which is simply proven by induction. 

Denote now the error in x by e = — — x, and let D = IJ — GA. Thus, 
from (3), 


De + Dette, 


h=1 


Restrict G by the condition that D? > 0 as p > ». Then, 


(4) e= 

h=1 

and this expression will be the basis for placing bounds on the error e. 

3. Requirements on D and First Uses of the Abmatrix. In Part 2 it was 
required that D? + 0 as p — ~. By means of the abmatrix, easy sufficient 
conditions for this convergence can be found. 

Let y = (g.), « = 1(1)”, be a column vector with g; > 0. It is always 
possible to find a scalar k, depending on , such that 


(5) a(D)y < ky. 
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It then follows that 
0 a(D*)y a(D?™)a(D)y ka(D?)y. 


Thus, by induction, 
€ a(D*)y kPy. 


Hence, if k <1, o~. 

In conclusion, (5) with k < 1 is a sufficient condition for the desired 
convergence. Clearly, the same arguments hold for a row vector 7’. 

Suppose then that k < 1 for some specified column vector +. Bounds 


on > D* may be found by column-wise examination. Consider > Dre; 
aban e; is the jth column of J. Then, 
h=1 h=1 
Since 
a(D)e; < fi, 
where f; is a scalar depending on 7, it follows that 


and for k < 1, : 
L h=1 i-k 
Similarly, for some specified row vector with k < 1, 
a[ < Gy. 
L h=1 


4. Error Bounds on the Approximate Solution. The application of (6) 
and (7) to stating error bounds on the approximate solution is immediate. 
It follows from (4) that 


(8) a) awe, 
h=1 
and (8) may be coupled with (6) or (7). It is true that (8) requires the matrix 
multiplication Dx, but considering that the larger matrix multiplication GA 
is required anyhow, it appears wise to benefit from the usually sharper error 
bounds given by (8) at the expense of some incremental work. In addition, 
(3) indicates that Dx is the first order correction to the approximate solution 
x, and thus Dx can also be used to improve x. 
If Dx is not computed, (4) may be used as follows: 


(9) al) <a[ ae 


and (9) may be coupled with (6) or (7). 
Finally, x may be taken to be a column of G, the approximate inverse, 
and thus, error in G may be also investigated. 
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5. Specific Example of Error Bounds on the Approximate Solution. As 
yet, no mention has been made of how y should be chosen once D is given. 
Generally speaking, the choice would be to minimize &, but this indeed is a 
serious problem.' Without looking for the optimum choice, inspection of D 


may give some ‘guide. However, a simple choice for y is e = z e;. For 
example, with D = (d,;), let 
k = max 9 
i i 
fi = max|dj;|. 
7 


It is then seen that 
e’a(D) < ke’, 


e’a(D) < fe’. 


fi 
«| i — Re? 


fi 


If k < 1, (7) gives 


whence from (9) 


(10) a(e,’e) 


e’a(x). 


Clearly, with the methods presented, many other error bounds such as 
(10) are possible. 
A. DE LA GARZA 
Carbide and Carbon Chemicals Company 
K-25 Plant, Oak Ridge, Tennessee 


1 The referee makes the following comments: 

“Tf all |d;;|>0, the optimum choice of y is unique and is the eigenvector ™; of a(D) 
whose components are all positive, and k is then the dominant eigenvalue \, of a(D). This 
follows from a lemma that, since all g; > 0, ; lies strictly between the minimum and maxi- 
mum of the ratios ¢;’a(D)/e;’y, unless the ratios are all equal (and hence equal to d;). The 
lemma is a slight extension of Theorem I of Haze. Perrect, ‘On matrices with positive 
elements,’ Quart. Jn. of Math., s. 2, v. 2, 1951, p. 286-290. 

“The vector a?(D)e, which is as mptotically a multiple of u; as p—> «©, may bea useful 
approximation to + for sufficiently p.” 
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1065[A].—H. S. Unter, ‘“‘Many-figure approximations for V2, ¥3, V4, and 
V9 with x? data,” Scripta Math., v. 18, 1952, p. 173-176. 


. Values of the constants mentioned in the title are given to 710, 709, 478, 
and 478 decimals respectively. 

The author reports errata in the values of the cube roots of 2 and 3 as 
given by Boorman.! These are correct only to 198D and 54D instead of the 
305D and 241D given. 

There is given also data on the distribution of the digits of each of the 
four cube roots. Values of x? and the associated probabilities are given for 
the first 50 k digits for k=1, 2, ---. No unusual distributions are revealed. 

D. H. L. 


1J. M. Boorman, “‘Square-root notes,” The Mathematical Magazine, v. 1, 1887, p. 208. 
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1066[D, R].—HeEtnz (1) Geoddtische Registertafeln. Schnell-Rechen- 
tafel. Finfstellige Winkelfunktionen. 400° mit direkter zweistelliger Inter- 
polation und Fehlergrenzen. Goslar-Harz, Selbstverlag, 1952, 60 p., 
16.2 X 24 cm., 9.80 Marks bound. (2) Geoddtische Registertafeln Schnell- 
Rechentafel. Sechsstellige Winkelfunktionen. 400° zur Maschinenrech- 
nungen. Goslar-Harz, Selbstverlag, 1953, 59 p., 16.7 K 25.5 cm., 14.00 
Marks bound. 


These very neatly prepared geodesy tables, each with many label flap 
guides, contain tables of sin, tan, cot, cos, for interval hundredth of a grade, 
differences being printed in red. In each is introductory material with 
illustrative examples. In (1) are a few useful supplementary tables; and in 
(2) are sine and cosine tables for [0(0*.1)100¢; 11D]. 

R. C. ARCHIBALD 


Brown University 
Providence, R. I. 


1067[F].—H. Gupta, “A table of values of N,(#),” East Panjab Univ., 
Research Bulletin, no. 20, 1952, p. 13-93. 


The table gives for each possible integer ¢ < 20000 all non-negative 
integral solutions x, y of 


(1) 


There is also given the sum function 
t 
N.(t) = 
j=0 


where m2(#) denotes the number of solutions (x,y) of (1), the solutions 
(xo, Yo) and (yo, xo) being counted as different when xo ¥ yo. 

This table is mentioned as existing in manuscript in MTAC, v. 3, p. 22, 
and was used to prepare the author’s table of N;(#), the number of solu- 


tions of 
Si. 
The author states that 
(2) N= 
0<i<4 


but this is in contradiction with the definition of N2(#) and its tabulated 
values. In fact the terms 
1+ [4] 


should be added to the right hand side of (2). 
The table can be used to find the exact number of lattice points inside 
or on a circle of radius r by means of the formula 


4N.(r*) — 4r — 3. 
Thus the number of lattice points inside and on a circle of radius 140 is 
4N,(19600) — 560 — 3 = 61529. 


The area of the circle is greater than this by the factor 1.000751. 
D. H. L. 
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1068[F].— Ove HEmMeER, On the Diophantine Equation y* — k = x*. Uppsala, 
1952, 101 p., 23.5 X 17.6 cm. 


The booklet is primarily concerned with the solutions of the equation 


y—-k=x* 
for positive k < 100. 

Table 1 (p. 94) purports to give all solutions for 47 of the 59 values of 
k < 100 for which solutions exist (10 of these were cases which were known 
previously), together with the number N of solutions. For k = 37, however, 
N is given as 2, but besides the two solutions which are listed, the solution 
(243, 3788) can be read from RoBINnson’s tables, which are on deposit in 
the UMT File [MTAC, v. 5, p. 162], and of which the author is apparently 
unaware. The discussion pertaining to the case k = 37 (p. 75-78) is therefore 
incomplete. In ten of the remaining 12 cases one solution is given and it is 
stated that N < 2. In the other 2 cases k = 63 and k = 76, two and one 
solutions are given respectively with the remark that “‘in all probability 
there are no more solutions.” 

There is only a one page discussion (p. 88-89) for negative k. A statement 
is made on p. 88 that y? + 56 = x* has no solutions. However, Robinson's 
table lists (18, 76) which is indeed a solution. Table 2 (p. 95) lists the 44 
values of k < 0 for which solutions exist, together with the known solutions 
and their number. In 22 of the cases all solutions are known, in the remaining 
cases one or two solutions are given. An extra solution can be read from 
Robinson’s table for — k = 7, 28, 39, 47, 53, 55, 60, and 63. 

Table 3 gives solutions of y? + 27k = x* for |k| < 50. These were ob- 
tained with the help of CassELs" table of rational solutions [MTAC, v. 4, 
p. 202]. All solutions are supposed to be given, except for 11 values of . 


For k = — 11, however, the number of solutions is given as 8, while a ninth 
solution (1362, 50265) can be read from Rostnson’s table. In two other 
cases, k = 24 and 37 which are not complete, Robinson gives an extra 
solution. 


Table 4 gives the discriminants —D, the fundamental rings and units e 
in the cubic field corresponding to 36 square free values of k between 2 and 
100. In 10 of these cases ¢ is not “definitely proved to be a fundamental unit.” 

Emma LEHMER 
16556 Chattanooga PI. 
Pacific Palisades, Calif. 


1J. W. S. CassEts, “‘ The rational solutions of the Diophantine equation y* = x* — D,” 
Acta Math., v. 82, 1950, p. 243-273. 


1069[F].—M. Krartcuik, ‘‘On the factorization of 2" + 1,” Scripta Math., 
v. 18, 1952, p. 39-52. 


This expository article contains tables giving the complete factorization 
of 2" + 1 for those values of m for which such information was known to the 
author. There are four tables. The first is for 2" — 1 and for 


nm = 1(2)99, 105, 107, 111(2)117, 123, 127, 135. 
The second is for 2" + 1 and for 
n = 1(2)99, 105, 111, 123, 135. 
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The third is for 2" + 1 and for 
n = 4(4)100, 108, 120. 
The fourth is for 2* + 1 and for 
n = 2(4)130, 138(4)162, 170(4)190, 198, 210, 222, 234, 250, 258, 270, 330. 


The tables separate the algebraic factors from the primitive ones. The 
fourth table has column headings for the two factors in 


These two factors are inadvertently interchanged for » = 4k + 2 = 94, 
114, and 150. 
An additional entry can be given, namely, 


2° + 1 = 3-179-6202089- 1858477404602617. 


The factor 6202089 was discovered by the SWAC. The large factor was 
proved prime by A. L. Brown and the reviewer independently. 
Other additional entries may be made from Notes 131, 138, 142, and 146. 
D. H. L. 


1070{I].—J. M. HAMMERSLEY, ‘“‘Lagrangian integration coefficients for dis- 
tance functions taken over right circular cylinders,” Jn. Math. Phys., 
v. 31, 1952, p. 139-150. 
The author considers the problem of evaluating the sixfold integral of 
the form 


I= f t(r)dv,dv2/r" 


where dv; and dv: are two volume elements, a distance r apart, in a cylinder 
of diameter a and length ca and ¢(r) is an arbitrary function of r. The 
integral is approximated by 


I = (a/2)*x*(cyi + ys) 
where 


f 


The auxiliary functions g, and k involve the hypergeometric function 
F(a, 1/2; 2; x) (a = — 1/2, —3/2). 
The integrals y; and y2 are to be evaluated by the Lagrangian integration 
formulas 
= 
ye = tat) H(é), 
the sums being taken over equally spaced points &. 
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The author tabulates G(¢), —H(é) and (for checking purposes) G(é) 
+ for 
= 0(.01)2(.05)3(.25)10; 8D 


with the last place uncertain. 
G. BLANCH 
NBSINA 


1071[K].—T. W. ANnpDErson & D. A. “Asymptctic theory of 
certain ‘goodness of fit’ criteria based on stochastic processes,’ Annals 
Math. Stat., v. 23, 1952, p. 193-212. 


The table contained in this paper (p. 203) gives the argument z to 5D 
which satisfies a;(z) = lim Pr (mw? = z) for a;(z) = .01(.01).99, .999. This 


is the limiting distribution for large samples of nw*, a statistic attributed to 
von Misss for testing the hypothesis that ” independent, identically dis- 
tributed random variables have a specified distribution function F(x). 
w* is defined by w? = /%. [F.(x) — F(x) PdF where F,(x) is the sample 
cumulative distribution function for a sample of size n. 
W. J. Drxon 
University of Oregon 

Eugene, Oregon 


1072(K].—Z. W. Birnsaum, “Numerical tabulation of the distribution of 
Kolmogorov’s statistic for finite sample size,’’ Amer. Stat. Assn., Jn., v. 
47, 1952, p. 425-441. 


Let Xi, Xe, ---, Xn be a random sample from a population having a 
continuous cumulative distribution function F(x). Define the empirical 
distribution function Fy(x) by Fyw(x) = j/N where j is the number of ob- 
served values of X1, ---, Xw which are less than or equal to x, and let 


Dy = l.u.b. | F(x) — Fr(x)|. 


Table I (p. 428-430) gives values to 5D of P(Dy < c/N) for N = 1(1)100 
and c = 1(1)15. Table II (p. 431) gives, for N = 2(1)5(5)30(10)100, values 
to 4D of the 95th and 99th percentiles of the distribution of Dy and com- 
pares these with the limiting values. 

FRANK MASSEY 
University of Oregon 

Eugene, Oregon 


1073[K].—A. H. Bowker & H. P. GoopE, Sampling Inspection by Variables. 
McGraw-Hill Book Co., Inc., New York, 1952. xi + 216 p., 15.6 X 23.3 
cm. Price $5.00. 


This book is built around a set of single and double variables sampling 
inspection plans of the non-central ¢ type. For normal variables the cases of 
known and unknown standard deviation, and two-sided as well as one-sided 
specifications, are included. No sequential variables plans are described, for 
reasons of incomplete theory and elaborate calculation. 
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For roughly corresponding discussion of single, double (and sequential) 
attributes sampling inspection plans, see Statistical Research Group, 
Columbia University, Sampling Inspection, New York, McGraw-Hill Book 
Company, Inc., 1948. 

The heart of this book lies in Chapter 11 (pages 125-141) on the mathe- 
matics of the plans, and in tables and charts (pages 164-209), together with 
a number of other tables and charts scattered throughout the text. In addi- 
tion, there is much material on sampling inspection generally, such as the 
installation and administration of sampling inspection plans, tightened and 
reduced plans, disposition of rejected product, non-normal populations, 
control charts in sampling inspection, comparisons of variables and at- 
tributes plans. 

A few of the principal tables: Table A (necessary for later tables): for 
each of fourteen inspection lot sizes and for each of three inspection levels, 
a sample size letter (42 entries, 14 different). Table B: for single sampling, 
for each of the fourteen sample size letters (each here corresponding to a 
single sample size) and for each of fifteen AOL classes, ranging unevenly 
from .024-.035 to 8.5—11.0, the value of k to 3D in the acceptance criterion 
=+ ks S Uorzé — ks = L, where U and L are the upper and lower specifi- 
cation limits respectively, and AQL, #, and s have the usual meanings. 
(The acceptable quality level, AQL, as here used, is the percentage of de- 
fective items in an inspection lot such that the inspection plan will accept 
95% of all lots submitted containing that percentage of defects.) Table C: 
as in Table B, for double sampling, showing ka, k-, and k;, all to 3D, corre- 
sponding to k in Table B, but for acceptance on the first sample, rejection 
on the first sample, acceptance on the second sample. Table D: for each of 
the fourteen sample size letters and for each of the fifteen AQL classes, the 
AQL and the LTPD. (The lot tolerance percent defective, LTPD, as here 
used is the percentage of defective items in an inspection lot such that the 
sampling plan will reject 90% of all lots submitted containing that percent- 
age of defects.) Table E (the major table of the book): for each of the four- 
teen sample size letters, the operating characteristic curves of each of from 
9 to 16 single and double sampling plans, each pair of plans showing the 
AQL class, the AOQL class, and the criteria for acceptance (163 curves in 
all). (The average outgoing quality limit, AOQL, is the maximum percent 
defective remaining in lots subjected to a plan if all rejected lots are screened, 
no matter what may be the quality of the lots submitted.) 

This work is a sturdy accomplishment by the authors and a credit to 
all who pioneered in this area. Among the latter are H. Romic, H. Dopce, 
N. L. Jounnson, B. L. WEtcu, W. J. JENNETT, C. C. Craic, W. A. WALLIs, 
K. J. ARNOLD, J. WoLFow!Tz, and the two authors. 

H. A. FREEMAN 


Massachusetts Institute of Technology 
Cambridge, Massachusetts 


1074[K].—K. FREUDENBERG, “Die Grenzen fiir die Anwendbarkeit des 
Gesetzes der kleinen Zahlen,’’ Metron, v. 16, 1951, p. 285-310. 


The principal objective of this paper is to determine the smallest sample 
size required in order to obtain a given degree of approximation to the 
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binomial distribution by means of the Poisson distribution and by means of 
the following, so to speak, second order approximating function 


_ (m — dj? — 


where s is sample size, \ = sq is the mean number of occurrences of an 
event of probability g in random samples of size s; \ is also the mean of the 
approximating Poisson distribution and m = an integer. The mean of By), m 
is also independent of s, and its variance is \(1 — As~). In spite of the fact 
that the zeroth moment of B,,,, is unity independent of s, it is not a distribu- 
tion, since some of its values can be negative. Since these negative prob- 
abilities are, however, exceedingly small, they do not vitiate the practical 
usefulness of B),, as an approximation to the Poisson and binomial distri- 
butions. In fact, the B,,,, require much smaller sample sizes for equivalent 
approximation to the binomial distribution than does the Poisson distribu- 
tion. The Poisson and binomial are in this notation denoted as follows: 


(2) = €"/m! 


an = - 


The following tables are of principal interest: on p. 307 values of sp 4 
and s,s as functions of \ = 3$(3)10, which are respectively lower bounds 
for s such that the deviations of A,» and By» from C\m =p = .001. s,4 
varies from 141.5 to 640.8, where correspondingly s,s goes from 7.8 to 78.8. 
On p. 308 s,,4 and s,,2 are given as functions of g = .01(.01).10, again for 
p = .001. Here s, 4 varies from 406.0 to 4912.2 while s,, 5 goes from 0 to 34.1. 

Auxiliary tables of some interest are found on p. 292-294. Aym, By,m, 
(Arxm — and (By m — Cym) are given for s = 50 and 100 and 
for\ = 3,m = 0(1)7, > 7;\ = 1, 3, m = 0(1)10, > 10;X = 2,m = 0(1)11, 
> 11; = 3, m = 0(1)12, > 12; = 4, m = 0(1)13, > 13. 

For the approximation B,,,, to be of practical use in sampling, extensive 
tables would be required; the present table has limited usefulness. The 
reader’s attention may be called to errors in formula (11) p. 289 and the 
subsequent argument. The right side of (11) should be \ and the subsequent 
argument can be corrected so as to yield this result. The findings of the 
paper are not affected by these errors. 

S. B. LITTAUER 
Columbia University 
New York City, New York 


1075[K].—C. H. GouLpEen, Meihods of Statistical Analysis. John Wiley & 
Sons, New York. 2nd ed., 1952, vii + 467 p., 15.2 XK 23.7 cm. $7.50. 


The tables in this volume present no novelty, consisting of areas and 
ordinates of the normal curve, percentage points for the Student-Fisher ¢, 
x?, and the variance ratio, values of the bias factor cz for the standard 
deviation in samples from normal, and an extract from the FisHER & YATES 
random numbers. 
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1076[K].—A. Hap, Statistical Tables and Formulas. London, Chapman & 
Hall; New York, John Wiley & Sons, 1952. 97 p., 27.8 X 21.9 cm., $2.50. 


The actual tables are preceded by 22 pages of pertinent formulas, ex- 
planations and illustrations. Tables I and II of ordinates and cumulative 
frequencies of the normal frequency function has the unusual feature of 
giving all values to 4S for the standardized argument u = 0(.01)4.99. 
Table III is a 3D probit table slightly abridged from Fisher & YarTeEs.' 
The usual percentage points of the Student-Fisher ¢ distribution are re- 
placed by “‘fractiles’” (Table IV) which are percentage points of the cumu- 
lative distribution. These are taken from the Fisher & Yates tables with the 
omission of 55% fractiles and the addition of 99.9% values as well as the 
set of values for 50, 80, 100, 200 and 500 degrees of freedom. Table V, frac- 
tiles of the x? distribution, is reproduced from the HALD & SINKBAEK table? 
and gives values of x¢? satisfying P(x? < x0?) = p to 1D for 100p = .05, 
1, .5, 1, 2.5, 5, 10(10)90, 95, 97.5, 99, 99.5, 99.9, 99.95 and the degrees of 
freedom f = 1(1)100. Table VI gives fractiles of x?/f to 4D for 100p = .05, 
A, .5, 1, 2.5, 5, 95, 97.5, 99, 99.9, 99.95 and f = 1(1)100(5)200(10)300(50)- 
1000(1000)5000, 10000. Table VII, fractiles of the variance ratio, x:*f2/x2*f1, 
is an extension both of the Fisher & Yates tables and of those of MERRING- 
TON & THompson.® Values are given to 3S for 100p = 50, 70, 90, 99.9, 99.95 
with f; = 1(1)10, 15, 20, 30, 50, 100, 200, 500, ©; f. = 1(1)20(2)30(10)60, 
80, 100, 200, 500, © and for 100p = 95, 97.5, 99, 99.5 with f; = 1(1)20(2)- 
30(5)50, 60, 80, 100, 200, 500, ©; f, = 1(1)30(2)50(5)70(10)100, 125, 150, 
200, 300, 500, 1000, «©. Table VIII for fractiles of the range in samples from 
a normal universe has been reproduced in part from Tippett‘ and PEarson.® 
Here values are given to 2D for 100p = .05, .1, .5, 2.5, 5, 10(10)90, 95, 97.5, 
99, 99.5, 99.9, 99.95 and the sample size m = 2(1)20. Values of the mean, 
standard deviation and their ratio for the range in samples from normal are 
given to 3D for the same set of sample sizes. 

The next two tables are devoted to functions useful in the estimation of 
the mean, £, and the standard deviation, ¢, of a one-sided truncated normal 
distribution with a known point of truncation according to methods given 
in HALp’s text® of which these tables form the companion volume. In the 
case only a sample of from the truncated range is used, for arguments y, 
calculable from the sample, z = f(y) gives an estimate of the standard- 
ized point of truncation, ¢; another function, s = Zg(z), where Z is the 
sample mean, estimates and estimates Table gives z to 3D for 
y = .55(.001).91, g(z) with its first differences to 4D for s = — 3(.1)2, and 
the variances and covariance of m'#/e and n's/o to 4S and their coefficient 
of correlation to 3D also for z = — 3(.1)2. For a one-sided censored normal 
distribution the sample consists of m — a known values falling to the right 
of the known truncation point and a of whose values no more is known than 
that they fall to the left of the point of truncation. In this case ¢ is estimated 
by z = f(h, y) where h = a/n and y is calculated from the known portion 
of the sample, ¢ is estimated from #g(k, z), where Z is the mean of the known 
portion of the sample and g(h,z) is simply obtained from the tabulated 
function y’(z) and h, and again # estimates ¢. In Table X, f(h, y) is given to 
3D for h = .05(.05).5 and y = .5(.005).6(.01).8(.05)1.5; y’(z) to 5S with 
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first and second differences, the same variances and covariances as before to 
4S and the same coefficient of correlation to 3D, all for z = — 3(.1)2. 

Table XI gives two-sided 95% and 99% confidence limits for the prob- 
ability 3 in the binomial distribution if the sample consists of x occurrences 
in ” trials. Values are given to 3D for x and m — x = 0(1)20(2)30(5)50, 60, 
80, 100, 200, 500, ~. Table XII gives 2 arc sin Vx to 4D for x = 0(.001)1. 
Tables XIII of log (m!) to 4D for m = 1(1)1000, XV of n? for n = 0(1)999, 
XVI of mt and (10m)! form = 1(.01)9.99, XVII of n— to 5D for m = 1(.01)9.99, 
and XVIII a 4D table of common logarithms are all readily found elsewhere. 


But Table XIV of log (2) to 4D for x = 1(1) (5) and m = 2(1)100 appears 


to be an extension of previous tables. And the final Table XIX is a new set 
of 15000 random digits compiled from drawings in the Danish State Interest 
Lottery. 
This is an excellent set of statistical tables containing a good deal of 
new material. 
ca 

1R. A. Fisher & F. YATEs, ae Tables for Biological, Agricultural and Medical 
Research. London & Edinburgh, 194 

2A, Hatp & S. A. SINKBAEK, table of percentage points ¢ x?-distribution,” 
Skandinavisk Aktuarietidskrift, v. 33, 1950, P.. 168-175. [MTAC, v. 

3M. MERRINGTON & C. M. THOMPSON, ‘ ‘Tables of percenta: tt. of the inverted 
beta A ) distribution,” Biometrika, v. 33, 1943, p. 73-88. [MTAC, v. 1, p. 78.] 

» TIPPETT, “On the extreme individuals and the range of samples taken from 

a norma population,” Biometrika, v. 17, 1925, p. 364-387. 

ARSON & H. O. HARTLEY, éThe probability integral of the range in samples 
of ’ obserration from a normal population, ” Biometrika, v. 32, 1942, p. 301-308. [M7 TAC, 
v. 
6 A. HALD, Statistical Theory with Engineering Applications. New York & London, 1952. 


1077(K].—T. Kiracawa, Tables of Poisson Distribution. Tokyo, Japan, 
Baifukan, 1951, 156 p., 18.7 X 26.0 cm., $3.50. 


The integer-valued random variable X is said to have a Poisson distribu- 
tion with parameter m if Pr (X = x) = e~™m*/x! = p(x, m). 

The tables give p(x, m) for x = 0(1)---, m = .001(.001)1 to 8D (Table 
I);m = 1.01(.01)5 to 8D (Table II); and m = 5.01(.01)10 to 7D (Table III). 
These tables were prepared directly from their definition. No interpola- 
tion methods are suggested. 

In the introductory material there is an English translation of a Japanese 
article by KiTaGawa on the subject of double sampling inspection plans 
which illustrates some of the uses of the above tables. 

In Table IV 95% and 99% confidence limits are given to 6D for the 
population parameter, m, with the sample mean # = .001, .005, .01, 
.02(.02).1, .2(.2)1(1)10 and the sample sizes 50(10)100(100)500, 1000, 
2000(2000)10000, 20000, 50000, 100000, when large sample theory is 
applicable. 

Kitagawa prepared his tables before it was possible for him to learn 
about Mottna’s tables.! In these (x, m) is given to 6D for m = .001(.001)- 


.01(.01).3(.1)15(1)100 together with P(c,m) = > p(x, m) to 6D for the 


same values of m and c = 0(1)153. Kitagawa’s volume does not contain a 
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table of either P(c, m) or 1 — P(c + 1, m). In about one hundred compari- 

sons between Molina’s and Kitagawa’s tables no discrepancies were found. 
I. R. SAVAGE 

National Bureau of Standards 

Washington, D. C. 


€. Mouina, Poisson’s Exponential Binomial Limit. New York, 1942. [MTAC, v. 1, 
p. 18. 


1078[K].—K. Kuntsawa, H. Maxase & H. Morimura, “Tables of confi- 
dence bands for the population distribution function. 1.,” Union of 
Japanese Scientists & Engineers, Reports of Statistical Application Re- 
search, v. 1, 1951, p. 23-44. 


Let (Xi, ---, X,) be a sample of independent variables with the same 
continuous cumulative distribution function F(x), and let N(z) be the 
number of X;’s which are =z. Setting F,*(x) = N(z)/n, the maximum D, 
of | F,*(z) — F(z)| is a random variable. Then Kotmocorov,! who intro- 
duced this statistic, showed that 1D, has the asymptotic cumulative prob- 


ability function L(z) = 1 — 2 © (— 1)*"e?#*#. This function was tabu- 
j=l 


lated by N. Smirnov? to 6D, for arguments z = .28(.01)3. In the present 
paper the authors publish a table of L(z) to 7D, for s = .250(.001)2.299, 
tabulating also the first differences. 
Z. W. BIRNBAUM 

University of Washington 
Seattle, Washington 

1A. Ko_mocorov, “Sulla determinazione empirica di una legge di distribuzione,”’ Inst. 
Ital. Attuari, Giorn., v. 4, 1933, p. 83-91. 

2N. Smirnov, “On the estimation of the discrepancy between empirical curves of 
distribution for two independent samples,” Bull. Math. de ? Univ. de Moscou, v. 2, 1939, 


~~ z= bs p. The table referred to was also reproduced in Annals Math. Stat., v. 19, 1948, 
p. 


1079[K].—F. J. MAssEy, JR., “Distribution table for the deviation between 
two sample cumulatives,’”’ Annals Math. Stat., v. 23, 1952, p. 435-441. 


Given two ordered samples of sizes N; and N, from cumulatives F;(x) 
and F;(y). Sy,(x) is the proportion of the members of the sample N; (¢ = 1, 2) 
less than x. Let d = max | Sy,(x) — Sw,(x)|. Table 1 (p. 436-439) presents 
all possible values of d, and a, for N;, N2 = 10, where 


a = Prob [d = Fi(x) = F;(x) }. 


The table is in terms of Ni, N2, 4, and a; h is an integer, such that d. = h/l, 
where / is the lowest common multiple of N; and Nz (kh SJ); and a is 
given to 5D. 

Table 2 (p. 440-441) is the same for additional selected values of N; 
and N, except that h usually is not carried out to /. 

R. L. ANDERSON 

University of North Carolina 
Raleigh, North Carolina 
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1080[K].—J. M. May, “Extended and corrected tables of the upper per- 
centage points of the ‘Studentized’ range,” Biometrika, v. 39, 1952, 
p. 192-193. 


This paper contains tables of the upper percentage points of the ‘Student- 
ized’ range g = (x, — %1)/s for the 5% and the 1% levels of significance for 
nm = 2(1)20 and »v = 1(1)20, 24, 30, 40, 60, 120, «© where s* is a mean square 
estimate of o” independent of the range, ” is the sample size for the range, 
and » is the number of degrees of freedom of s. Most entries are given to 2D. 
None contains fewer than 3S. These tables were computed by numerical 
quadrature as a supplement to tables computed earlier by PEARSON & 
Hart Ley! with the aid of approximation formulas which were found un- 
reliable for small » and large g. The Pearson-Hartley tables did not include 
entries for vy < 10 and the authors pointed out that entries exceeding 6 
should not be considered reliable for y = 10. For additional tables of the 
‘Studentized’ range see the review of PILLat’s paper in this issue [RMT 1083]. 
A. C. CoHEN, Jr. 
University of Georgia 

Athens, Georgia 


1E. S. Pearson & H. O. HaRTLEy, _ of the probability integral of the ‘Student- 
ized’ range,” Biometrika, v. 33, 1943, p. 89- 


1081[K].—J. Mosman, “Testing a straggler mean in a two-way classifica- 
tion using the range,’’ Annals Math. Stat., v. 23, 1952, p. 126-132. 


In 1948 Narr? obtained the distribution of the Studentized extreme devi- 
ate from the sample mean, i.e., he obtained the distribution of the sample 
statistic 

U, = (x, — Z)/s, or u,’ = ( — x)/s, 


k 
where x} S Sx, =k“ x; and s, is an independent estimate 


i=1 

of the standard deviation based on v degrees of freedom. In a two-way 
classification in the analysis of variance, therefore, s, can be taken as the 
residual standard deviation and the above test can be used to judge whether 
a row or column mean is high or low relative to the grand mean. In 1950 
PATNAIK? published a paper which indicated that the distribution of the 
mean range can be approximated with sufficient practical accuracy by a x? 
distribution, provided an appropriate scale factor is used and an equivalent 
(usually fractional) number of degrees of freedom employed in the conver- 
sion. Tables were published in Patnaik’s article which give both the appro- 
priate scale factor and the equivalent number of degrees of freedom for the 
mean range. These two results have apparently led the author of the present 
paper to recommend a test for detecting a “‘straggler”’ row or column mean 
in a two-way classification using in effect the mean range instead of the 
standard deviation in Narr’s test. As a matter of fact, the test proposed by 
MosHMAN is based on the sample statistic 


r r ik k 
g= i=1 i=1j=1 | j=1 f=1j=1 
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The numerator of the g statistic hence involves only row or column totals, 
the grand total, the number of rows and columns and a function W,,, of 
certain “ranges.’’ This latter factor is in fact the sum of the “column” 
(‘‘row’’) ranges of the “row” (“‘column”’) residuals. Appropriate tables for 
the 1% and 5% significance levels of the statistic, g, are given to 3S in the 
article for all cases where the number of rows or columns does not exceed 
nine. A quick answer concerning whether a particular row or column mean 
in a two-way classification is a straggler can therefore be obtained directly 
from these tables. For another method see TuKEyY.* 
F. E. Grusss 

Ballistic Research Laboratories 
Aberdeen Proving Grounds, Maryland 

1K. R. Narr, “The distribution of the extreme deviate from the sample mean and its 
Studentized form,” Biometrika, v. 35, 1948, p. 118-144. 

2 P. B. Patnatk, “The use of the mean range as an estimator of variance in statistical 
tests,’ Biometrika, v. 37, 1950, p. 78-87. 

3 J. W. Tuxey, “Comparing individual means in the analysis of variance,” Biometrics, 
v. 5, 1949, p. 99-114. 


1082[K].—D. N. Nanna, “Probability distribution tables of the larger root 
of a determinantal equation with two roots,” Indian Soc. Agr. Stat., Jn., 
v. 3, 1951, p. 175-177. 


Define the distribution function P(x) by 


2P(x)I(1, 2m + 2, 2n + 2) = 21 (x, 2m + 2, 2n + 2) 


— — (x, m + 1,” + 1) 
where 


The author tabulates and P-*(.99) to 2D for m = 0(})2,” = 4$(4)10. 
These tables provide tests for a number of statistical hypotheses, includ- 
ing that of the equality of two covariance matrices of bivariate normal 
populations. 

J. L. Hopcgs, Jr. 
University of California 
Berkeley, California 


1083[K].—K. C. S. Prixat, “On the distribution of the ‘Studentized’ range,” 
Biometrika, v. 39, 1952, p. 194-195. 


This paper contains tables of the lower percentage points of the ‘Student- 
ized’ range g = (x, — x1)/s at the 5% and the 1% levels of significance for 
n = 2(1)8 and » = 1(1)9, where the notation is that of RMT 1080. All 
entries are given to 2D. These results further supplement those of PEARSON 
and Hart ey.! An editorial note explains that PILLAI also computed certain 
upper percentage points which were not published separately since his paper 
was submitted while May’s (see RMT 1080) more extensive program of 
computation was in progress. Pillai’s upper 5% and 1% points computed 
for 1 Sv» =4 and 2 =” S 8 are included in the tables by May, and ac- 
cording to May, agreed with his calculations to within a unit of the last 
figure quoted by Pillai. In order to perform the necessary calculations, the 
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author determined the distribution of g explicitly in the form of an in- 
finite series. 

A. C. COHEN, JR. 
University of Georgia 
Athens, Georgia 


1E. S. Pearson & H. O. HartTLeEy, yw of the probability integral of the ‘Student- 
ized’ range,” Biometrika, v. 33, 1943, p. 89- 


1084[K].—ALan Stuart, “The power of two difference-sign tests,’’ Am. 
Stat. Assn., Ju., v. 47, 1952, p. 416-424. 


The author considers a distribution-free test of serial independence of NV 
(unequal) observations with its N — 1 first differences ordered in time, a 
test proposed by Moore & WALLIs.! 

w-1 


Define D = > D; where D; is unity when the ith difference is positive, 


i=1 
and zero otherwise. The null hypothesis is that the observations come from 
the same continuous population, and under this hypothesis the mean and 
variance of D:is derived and its asymptotic normality for large N proved. 
The alternative hyputhesis, assuming the observations equally spaced, is 
Ye = a + BR + x, where x; is a normal variate with mean 6 and variance o?. 
The mean, variance, and the asymptotic normality under the alternative 
hypothesis is derived. Assuming exact normality of the statistic D under 
both the null and alternative hypotheses, the author tabulates the power 
of the D-test to 2D at the 95% level of significance for N = 15, 25, 50, 75, 
100 and R = B/oV2 = .1(.1).5, 1. He concludes for N = 25, R =.5, and 
N 2 75, R = .3, that the loss of power as compared with the best available 
test is at most 8 percent. He next considers the case of two series of NV un- 
equal observations in time, and to investigate the correlation between them 
defines D; for the first series as before, and similarly ¢ d; for the second series. 


Let C; = D,d;, and the correlation statistic C = = C;. Under the null 


hypothesis that the N pairs of observations are drawn from the same con- 
tinuous bivariate population in which the two variates are independent he 
finds the mean variance of C and proves its asymptotic normality. The 
alternative hypothesis is that the N pairs of observations are drawn from a 
normal bivariate population with p ¥ 0, and under this hypothesis the mean, 
variance of C and its asymptotic normality is derived. The author provides 
a table for the power of the C test to 2D against this alternative (under the 
assumption again of exact normality) at the 95% level for p = .1(.1).7 
and N = 50, 100, 200, and 400. This is compared to the power of the usual 


Z test, Z = } log : 3 4) at the 95% level for p = .1(.1).5, NW = 50, 100, 


200, 400, with results given to 2D. For p = .5, N = 100, p = .4, N = 200; 
or p = .3, N = 400, he concludes that the use of the C test is satisfactory 
since at most 8% of the power is lost. 


L. A. AROIAN 
Hughes Aircraft Company 
Culver City, California 


1G. H. Moore & W. A. Wattis, “Time series significance tests based on signs of 
differences,”’ Amer. Stat. Assn., Jn., v. 38, 1943, p. 153-164. 
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1085[K].—T. J. Terpstra, ‘‘A confidence interval for the probability that 
a normally distributed variable exceeds a given value, based on the 
mean and the mean range of a number of samples,” Applied Scientific 
Research, A, v. 3, 1952, p. 297-307. 


The author employs PatNarx’s' method to find the confidence intervals 
described in the title. For / samples of size k each the author gives four charts 
from which one can read the central 90% confidence limits on the vertical 
axis with the ratio of the numerical difference of the given value from the 
mean of the total sample to the average range on the horizontal axis, for 
1 = 4(2)10 and k = 1(1)5. 


1P. B. Patnatk, “The use of the mean range as an estimator of variance in statistical 
tests,” Biometrika, v. 37, 1950, p. 78-87. 


1086[K].—J. WisHart, ‘Moment coefficients of the k-statistics in samples 
from a finite population,” Biometrika, v. 39, 1952, p. 1-13. 


TuKEy! in order to simplify the problem of finding moment coefficients 
of Fisher’s k-statistics in samples from a finite population introduced the 
sample statistic K,,... which has the property that its expected value is the 
same function of the N members of the population that &,,... is of the ” 
members of the sample. The author first gives a conversion table for finding 
the k,,... statistics in terms of the augmented monomial symmetric functions 
of Davin & KENDALL? and vice versa through order 6. Then to facilitate the 
derivation of further sampling formulae he gives a complete table of products 
and powers of k-statistics in terms of the ,,... through order 6. For orders 
7 and 8 he gives the 14 and 21 basic formulae respectively, the remaining 
results of those weights then being derivable by algebraic substitution. 
Illustrations of the application of these tables are given. 

oni ‘St W. TuKey, “Some sampling simplified,” Amer. Stat. Assn., Jn., v. 45, 1950, p. 


2F. N. Davip & M. G. KENDALL, “Tables of symmetric functions—part I,” Biometrika, 
v. 36, 1949, p. 431-449. [MTAC, v. 4, p. 146.] 


1087[L].—BritisH ASSOCIATION FOR THE ADVANCEMENT OF SCIENCE, Com- 
mittee on Mathematical Tables, Table of Bessel Function, Part I1. Cam- 
bridge, 1952. xl + 255 p., 21.5 X 28 cm. 


This is the long-awaited ‘‘second volume.”’ The functions tabulated are 
Jn(x), Ya(x), In(x) = i-*J,(ix) and K,(x) = (ix) + ¢Y,(éx)] for 
n = 2(1)20 (the first volume! gave values of the functions for = 0, 1) and 
for x = 0(.1 or .01)10(.1)20 with 8? or 6,,?. Tables for J, are'8D whereas 
those for Y,, J,, and K, are to 8 figures. In regions where Y, J, or K is 
not convenient for interpolation, the functions x*Y,, x*K,, eI, or eK, 
are tabulated. Finally J,, Y., In, and K, are tabulated for » = 0(1)20, 
x = 0(.1)25, to 10D for J,, 10 figures for Y,, I,, Kn. 

The tables are, in part, based on those of ALprs? but are very considerably 
expanded and reworked. The checking has been carried out with the 
scrupulousness typical of this series of tables. 
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Bessel functions are like the weather, they nearly always intrude some- 
where in any research in theoretical physics or applied mathematics. Like 
the weather also, for a long time very little was done about assembling a 
set of tables of the functions complete enough to satisfy most needs. Few 
of the calculations of acoustical and electromagnetic wave radiation, propa- 
gation or scattering can be made with the sole use of Bessel functions of the 
first kind, of real argument, so attempts to use the remarkable HARVARD 
CoMPUTATION LABORATORY tables are usually reminiscent of attempts to 
skate with only one ice-skate. Here, in the present BAASMTC volume, we 
have both solutions tabulated, for both real and imaginary values of the 
argument, with enough significant figures and enough range of the argument 
and order so that most calculations can be carried through. The Bessel 
function Panel responsible for its production—W. G. BickLey, L. J. 
Comrie, J. C. P. D. H. and A. J. THompson—are to be 
congratulated. Incidentally, this is the last volume of the BAAS tables; 
the Committee has been reconstituted as the Royal Society Mathematical 
Tables Committee. 

As far as the needs of the theoretical physicist and engineer are con- 
cerned, this volume and the earlier one! will suffice, in the majority of cases. 
In a few scattering calculations the functions J and Y for m = 20(1)50 for 
x = 0(.1)25 would be useful, for in scattering calculations, orders up to 
about twice the argument, still contribute in the result. In some electro- 
magnetic and acoustical problems, functions of a complex argument are 
required. These are available in the NYMTP tables‘ for » = 0 and 1. It 
would be useful to have similar tables for m = 2(1)10, but such needs are 
not urgent. Tables available and desired for half-integer orders will be dis- 
cussed in RMT 1092. 

The volume under review here is clearly printed and each page is clearly 
marked for ready reference. It is a pleasure to use the volume. 


Puitie M. Morse 
Massachusetts Institute of Technology 
Cambridge 39, Mass. 


1BAASMTC, Bessel Functions, Part I. Cambridge, 1937, xx + 288 p. 
2W.S. Atpis, Roy. Soc. Proc., v. 64, p. 203, 1899, and v. 66, p. 32, 1900. 
8 HARVARD COMPUTATION LABORATORY, Tables of Bessel Functions of the First Kind for 
n = 0 to 135, vols. III to XIV, 1947-1951, Harvard U. Press, Cambridge, Mass. 
‘ NYMTP, Tables of Jo, Ji and Yo, Vi for Complex Arguments, 1943 and 1950, Columbia 
University Press, New York. 


1088[L].—H. Koser, Dictionary of Conformal Representations. New York, 
Dover, 1952. xvi + 208 p., 15.2 X 23.8 cm. Price $3.95. 


The book represents a formidable achievement in bringing together a 
wide variety of explicit conformal maps useful to physicists, engineers and 
computers. Particularly notable are the numerous clear diagrams, which 
allow the reader to locate quickly a required mapping. Most of the trans- 
formations are by elementary functions, but short theoretical discussions of 
the Schwarz-Christoffel transformation and of the Schwarz principle of 
reflection indicate in part the derivation of more complicated formulas. 
Proofs are altogether omitted due to the dictionary character of the presen- 
tation. Mappings are classified according to the nature of the analytic func- 
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tions involved, with chapters on bilinear transformations, algebraic func- 
tions, trigonometric and exponential functions, the Schwarz-Christoffel 
transformation, and elliptic and modular functions. Although a geometrical 
index is listed as topological, almost no reference is made to Riemann sur- 
faces. This is particularly to be regretted, since it is quite common in 
engineering problems, such as those involving the hodograph method, to 
ignore the usefulness and necessity of exploiting simple conformal mappings 
onto Riemann surfaces. Accuracy and a notable absence of misprints, of 
which the reviewer found only one, should make the handbook of permanent 
and growing value to the applied scientist. 

P. R. GARABEDIAN 
Stanford University 
California 


[Editorial Note. A preliminary version of this work was reviewed in MTAC, v. 2, p. 
296-297, v. 3, p. 103, 368]. 


1089[L].—N. H. Marcu, ‘‘Thomas-Fermi fields for molecules with tetra- 
hedral and octahedral symmetry,’’ Cambridge Phil. Soc., Proc., v. 48, 
1952, p. 665-682. 


Two types of solutions of the differential equation 
Pp 
— = x—ig! 
da * 


are tabulated in this paper. One of these is of the form 
= 1+ + ax! + ag? + --- 


and is fully determined once az is given. Table 2 (p. 670) gives as, ---, @11 
in terms of a2, and Table 1 (p. 668) values of ¢ for — a2 = 1.55, 1.56, 1.567, 
1.5748, 1.581, 1.5849, 1.5853, 1.5866, 1.58756, 1.58764. x varies from 0 to 
10.74. The solutions were started by using the above series and then con- 
tinued by numerical integration of the differential equation. Table 1 also 
gives in each case the point x» the tangent at which passes through the origin, 
with the corresponding ¢(x9). It is suggested that a classification of the 
solutions by $(xo) is more reliable than that based on az. 
The second type of solutions is of the form 


144 #F; 
1— +2-24...). 


Table 3 (p. 670) gives expressions of F2, ---, Fio in terms of F;. It turns out 
that all solutions of this type may be derived from two “master solutions” 
in which F, is, respectively, 1 or — 1. These two niaster solutions and their 
first derivatives are tabulated in Tables 4 and 5 (p. 671, 672) for values of 
x which range from .1375 to 20.2 in Table 4 (F; = 1), and from .5375 to 
19.4 in Table 5 (F; = — 1). 

In the paper there are also references to solutions computed previously.! 

A. E. 


1R, P. Feynman, N. MEtropouis, & E. TELLER, “Equations of state of elements based 

on the eee Fermi-Thomas Theory, ” Phys. Rev., s. 2, v. 75, 1949, p. 1561. 
J. C. Stater & H. M. Krutter, “The Thomas-Fermi method for metals, ” Phys. Rev., 
s. 2, v. 47, 1935, p. 559. 


me- 
ike 
ga 
‘ew 
pa- 
the 
to 

we 
the 
ent 
ssel 

be 
ical 
ses. 
for 
are 
are 
dlis- 
l for 
abia 
ra 
and 
lich 
ins- 
s of 
of 
las. 
en- 
inc- 


100 RECENT MATHEMATICAL TABLES 


1090[L, M].—W. Meyer zur CAPELLEN, Integraliafeln. Sammlung un- 
bestimmter Integrale elementarer Funktionen. Springer-Verlag. Berlin/ 
Géttingen/Heidelberg, 1950. viii + 292 p., 23.8 X 16.5 cm. 


This book contains probably the largest published collection of indefinite 
integrals whose integrands are combinations of elementary functions. It 
was undertaken in 1942 and completed in 1944, publication being delayed 
by difficulties of the post-war period. There are about 3000 entries, most of 
them giving the final answer, and some of them stating recurrence relations 
or other auxiliary formulae for the integrals in question. Many special cases 
of more general results were included, thus enabling the user to get the final 
answer quickly. For instance, in section 2.2.2, there are recurrence relations 
for the integral 


f x?(a + bx)"/™dx m,n, p integers, 


and there is also an explicit formula for positive integer p. Nevertheless, 
there are about 7 pages of some eighty particular cases listed, for special 
values of m, n, p—a great convenience for the user. 

Contents. Chap. 1. Preliminaries (5 p.). Rules for differentiation and 
integration. Chap. 2. Integrals of algebraic functions (138 p.). 2.1. Rational 
integrands. 2.2. Irrational integrands leading to elementary functions. 
2.3. Algebraic integrands leading to elliptic integrals. Chap. 3. Integrals of 
transcendental functions (80 p.). 3.1. Exponential and logarithmic functions. 
3.2. Trigonometric and inverse trigonometric functions. 3.3. Hyperbolic 
and inverse hyperbolic functions. Chap. 4. Products of algebraic and trans- 
cendental functions (45 p.). Arranged according to the transcendental part, 
similar subdivision to that of Chapter 3. Chap. 5. Products of transcendental 
functions (11 p.). 5.1. f-g(x) In xdx. 5.2. fe*g(x)dx. Addenda (3 p.). Chap. 6. 
Supplementary material (9 p.). Constants, expansions, definitions of some 
higher transcendental functions. Chap. 7. References (19 items, 1 p.). 

The classification and arrangement of such a large number of integrals 
presents very formidable difficulties, and only continued and frequent use 
of the tables will show the excellence, or the imperfections, of the scheme 
adopted here. The author’s aim was to provide tables which can be used 
with ease by a non-mathematician, and the first impression is decidedly 
favorable. There is a very detailed table of contents which enables one to 
locate any integral printed in the book within four or five pages. One will, 
of course, encounter integrals which are not in the book, yet can be trans- 
formed into integrals printed in the book; but a generous duplication of 
several forms of the same integral will greatly reduce the number of these 
occurrences. 

The only drawback the reviewer found in using these tables is the absence 
of a list of notations and abbreviations. The user not acquainted with 
German books will find out sooner or later that Sin means sinh and ArSin 
means sinh (although these symbols are defined on p. 286, 287 under hyper- 
bolic functions, meeting them on p. 77 as integrals of irrational functions 
gives no clue as to the place where one should look for their definition). 
Ad hoc notations are sometimes explained where they occur (for instance, 
x on p. 63 in 2.3.1.2.1.3.1) or in a footnote on the same page on which they 
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are used (as U, V, W, on p. 68); yet sometimes one has to turn several pages 
until one finds an explanation (the same V, W, are used on p. 76, and one 
has to go back to p. 68 to find a footnote referring one to nos. 2.2.1.1.1 and 
2.2.1.1.2 which are then found on p. 73, 74). Also, the numbering of the 
entries (samples above) is cumbersome, if systematic, and makes it incon- 
venient to give specific references. 

In all-other respects this is an eminently useful book, and the printing 
is excellent. 

A. E. 
[For errata in this work see in this issue MTE 223. ] 


1091[L].—NBS Applied Mathematics Series, No. 17. Tables of Coulomb wave 
functions. U. S. Government Printing Office, Washington, D. C., 1952. 
xxvii + 141 p., 20.5 X 27 cm. Price, $2.00. 


Coulomb wave functions are confluent hypergeometric functions, but 
the normalisation and notation usual in quantum mechanics differ from the 
standard notations in mathematical texts. 

The “regular Coulomb wave function,” Fz(y, p) is that solution of the 
differential equation 


@y 


which is regular at p = 0 and is so normalised that 
Fr(n, p) ~ sin (p — 9 In 2p — + ox) 
as p— ©, where oz = arg I'(L + 1 + im). The notation 
is frequently used, where 
Co(n) = — L(2L + 1)Ci(m) = (L* + 


and ;(n, p) is a power series in p. 

Coulomb wave functions have been computed by several authors: on p. 
152 of the recent tables by BLocu et al (MTAC, v. 6, p. 92), there\are references 
to available tables. The volume under review contains the most extensive, 
and the most systematic, tables, and it largely (but not entirely) supersedes 
the earlier tabulations. 

Table I (p. 1-111) gives 7D values of @z(n, p), and of as many reduced 
derivatives 


dy $1(n, p) 


as are needed for interpolation. In this table, y = — 5(1)5, L = 0(1)5, 10, 
11, 20, 21, and p = 0(.2)5. 

Table II (p. 113-127) gives 10D values of the real part of I’(1 + in)/ 
T'(1 + in), with modified second central differences, for » = 0(.005)2(.01) 
6(.02)10(.1)20(.2)60(.5)110. 


/ = 
| 
’ 


102 RECENT MATHEMATICAL TABLES 


Table III (p. 129-135) gives 8D values of oo, with modified second cen- 
tral differences, for 7 = 0(.01)1(.02)3(.05)10(.2)20(.4)30(.5)85. 

Table IV (p. 137-141) gives 8S values of Co() for 7 = 0(.01)5(.05)10. 

The Foreword contains an account of the applications of Coulomb wave 
functions, and the Introduction (by M. ABRAMOWITZ) gives definitions of, 
and formulas, Bessel function expansions, and asymptotic expansions for, 
Coulomb wave functions. The Introduction also describes the method of 
computation, and gives instructions for interpolation and for the use of 
recurrence relations. 

In view of the fact that this volume deserves, and is likely, to become 
the principal work of reference for Coulomb wave functions, without en- 
tirely superseding some of the earlier tables, it is regrettable that there 
is no systematic account of all (although there are references to some) 
available tables. 

A. E. 


1092[L].— Roya. Society MATHEMATICAL TABLES COMMITTEE, Short Table 
of Bessel Functions Ins3, Kn44. Cambridge, 1952, 20 p., 21.5 X 28 cm. 


This is a preliminary table of spherical Bessel functions of imaginary 
argument. The functions actually tabulated are 


x4 3 (x) for n = 0(1)10 and x = 0(.1)5, 

(2/m)x"*4Kn44(x) for n = 0(1)10 and x = 0(.1)5, 
€*1n44(x) for n = 0(1)10 and x = 5(.1)10, 
(2/m)e*Kn4s(x) for n = 0(1)10 and x = 5(.1)10, 


in most cases to eight significant figures, together with 6,,?, for interpolation 
by Everett’s formula, when J,(x) = 7-*J,(ix) and 


These tables begin the completion of the spherical Bessel functions as 
the table in RMT 1087 completes the main tabulation for cylindrical 
Bessel functions. The NBS tables! have covered the two independent solu- 
tions, for real values of the argument, for m = 0(1)30 and x = 0(.01 or .1)25 
to 8 or 10 significant figures. The present table is a start toward a similar 
table for imaginary arguments. 

Spherical Bessel functions of imaginary argument are useful in the 
quantum mechanical calculation of atomic and molecular properties. In 
particular, an atomic wave function for a bound state, about one center, 
when expanded about another center, appears as a sum of such Bessel func- 
tions, multiplied by spherical harmonics.? Calculations of interatomic forces 
in a molecule have been quite difficult to carry out, chiefly because of the 
lack of adequate tables. Now, with the advent of high speed computing 
machines, to remove the burden of routine details, such computations 
should become relatively easy to perform and our insight into molecular 
structure can begin to be satisfactorily complete. Values of J,4;(x) and 
Kaii(x) for » ~ 0(1)20 and x ~ 0(.1)20 will be needed as a basis for the 
calculation of the component integrals. The present table is a promising 
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start in this direction; let us hope that the RSMTC, or some other group, 
will carry on. 
M. Morse 
Massachusetts Institute of Technology 
Cambridge 39, Mass. 
 NBSCL, Tables of Spherical Bessel Functions. Vols. | and II, 1947, Columbia University 
Press, New York. 


3 See, for example, N. Rosen, “Interaction between Atoms with s-Electrons,” Phys. 
Rev., v. 38, p. 255, 1931. 


1093[V].— AERONAUTICAL RESEARCH COUNCIL, A Selection of Tables for Use 
in Calculations of Compressible Airflow. Prepared by the Compressible 
Flow Tables Panel (L. RosENHEAD, Chairman; W. G. BickKLey, C. W. 
Jones, L. F. Nicnotson, C. K. THORNHILL, R. C. ToMLinson) of the 
ARC. Oxford, 1952, viii + 143 p., 21.5 X 27.0 cm. Price $8.00. 


The earlier and more familiar tables for compressible flow in this country 
are H. W. Emmons’ “Gas Dynamics Tables for Air’’ (Dover publications, 
1947, MTAC, v. 3, p. 36) and J. H. Keenan & J. Kaye’s “Gas Tables” 
(Wiley, 1945, MTAC, v. 2, p. 20). Compared with these tables, the present 
publication is certainly more detailed and more accurate. It is stated that 
the two earlier tables have numerous errors, mostly affecting the final place 
of decimals. As in Emmons’ table, the value of the ratio of specific heats, +, 
is chosen to be 1.4. The authors try to remedy this limitation in part by 
listing the derivatives of a few important quantities with respect to y at 
y = 1.4. But nevertheless, this restriction will make the table unsuitable 
for calculations of gas flow at higher temperatures, such as the flow in high 
pressure compressors and gas turbines. 

There are three main groups of tables: The first group is concerned with 
isentropic one-dimensional flow. The argument used is the local Mach 
number M, the ratio S of speed to the critical sonic speed, or the square of 
the ratio of speed to maximum speed r. The M-tables have the ranges 
M = 0(.01).5(.001)1; 1(.01)5. The S-tables have the ranges S = 0(.01)2.44 
and ¥6; 1(.005)1.5/.01)2.44 and V6. The r-table has the range r = 0(.01)1. 

The second group of tables is for the solution by the characteristics 
method. The argument used is either the angular characteristics variable w 
or the Mach angle yu. The ranges of these tables are w = 0°(.5°)100°(1°)130°; 
0°(.001°).05°; 0°(.01°)1°, 0°(.1°)30°. » = 0°(1°)90°. 

The third main group of tables gives the ratio of quantities after the shock 
to that before the shock. The argument of the tables is either the Mach 
number M, or the speed ratio S; before the shock. For normal shock waves, 

= 1(.01)5; S; = 1(.01)2.44 and V6. For oblique shock waves only the 
maximum flow deflection angle and angle for sonic flow after shock are given 
for M, = 1(.01)2(.1)5(1)20. 

Other tables in this volume are tables for obtaining Mach number from 
the dynamic pressure ratio or the Pitot pressure ratio, tables of pertinent 
powers of x and (1 — x*), tables for calculating the Reynolds number, the 
pressure coefficient, etc. 


H. S. Tsien 


California Institute of Technology 
Pasadena, California 
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In this issue references to errata have been made in RMT’s 1065, 1068 
and 1069. 


221.—L. A. Livsrernix, I. fa. Axususxit, & V. A. Dirxin, Tablifsy Besse- 
levykh Funkisit. (Matematicheskie Tabliisy, no. 1.) [Tables of Bessel 
Functions (Mathematical Tables, v. 1) ] 1949. [MTAC, v. 5, p. 25.] 


In addition to the errors noted in MTAC, v. 5, p. 25, the following errata 
were discovered by inspection and checking of differences. 


x Jo(x) Ji(x) 
; for read for read 
0.167 9920399 9930399 
0.860 3914929 3914529 
1.678 694 674 
1.944 2592767 2562767 
2.171 1265437 1265435 
2.343 5119 5319 
2.983 3457353 3453753 
3.345 2000896 2020896 
3.662 608 698 
3.939 425404 425434 
3.960 506903 506953 
4.099 8887728 3887728 
4.316 4769 1769 
4.326 0801626 1801626 
6.755 791 789 
7.046 —0.0091064 -+0.0091064 
7.750 1916006 1916026 
8.141 2501847 2521847 
8.674 2707632 2707622 
9.849 2247 2447 
9.896 696 694 
10.028 2606 2506 
11.179 2914 2014 
12.350 1356 1354 
12.930 2040869 2014869 
13.152 2151370 2151350 
13.340 25587 35587 
13.575 2115728 2115428 
13.908 1179283 1179223 


14.317 1210012 1215012 
14.557 760050 764050 


14.562 750175 754175 
14.577 720459 724459 
14.837 190310 194310 


15.047 230301 238301 


104 
p. 
11 
25 
46 
51 
54 
67 
2 81 
87 
89 
94 
94 
143 
148 
162 
170 
181 
204 
205 
208 
231 
> 254 
266 
271 
274 
279 
286 
299 
299 
299 
304 
308 
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x Jo(x) Ji(x) 

p. for read for read 
310 15.128 1998981 1998921 
311 15.173 1972955 1972925 
328 16.260 1920047 1921047 
332- 18.490 162 164 
345 25.000 101 99 


Three dozen values in which the errors amounted to one unit in the last 
place have not been included in this list. 
S. A. JOFFE 


515 W. 110 Street 
New York 25, N. Y. 


222.—N. W. McLacuaian, P. HumBert, & L. Pot, Supplément au Formu- 
laire pour le Calcul Symbolique. Mémorial des Sciences Mathématiques, 
fasc. 113, 1950. 


In addition to the errata published in MTAC, v. 6, p. 101, the following 
errors occur. 


. 5, item 2; for Ei(— x) on the r.h.s. read Ei(x). 

. 5, line 10; imsert ) after x. 

. 6, formula 7; for the last + read —. 

. 7, formula 3; for x* read x~*. 

10, second item for p. 12; for o(log p) read o(log p) — f(0). 

10, second item for p. 14; for — 1 read — p. 

. 11, item for p. 34; read delete formulas 3-10. 

. 11, item for p. 47; for — a*/s* read — a?/s*. 

. 11, item for p. 48; for I,(p?/2) read I,,(2/p*) as in first edition. 

. 14, formula 10; for p* read p. 

. 15, formula 8; lower limit in second integral should be p. 

. 16, formula 8; for 2" read 2/2, and for x/2V¢t read x/ V2t. 

. 22, sixth formula from foot; in denominator for p read p’. 

. 26, formulas 4, 5; delete. 

. 27, last formula; replace |.h.s. by log {(¢ + 1)/(¢ — 1)}, and r.h.s. by 
e-?(log 2p + y) — e? Ei (— 2p),¢ > 1. 

28, formula 1; multiply the r.h.s. by 2. 

. 31, formulas 3, 8; the centre sign on the I.h.s. should be +. 

. 43, formula 9; on the r.h.s. for (ab)?™+"(2b)**+! read 2#tq2etepiettDtm, 

. 43, in the last two integrals, for dx read dt. 


I am indebted to A. ErpELY1 for some of these corrections. 
N. W. 


Vizard & Co. 
51 Lincoln’s Inn Fields 
London W.C.2, England 
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223.—W. MEYER zUR CAPELLEN, Integraltafeln, Springer, Berlin, 1950 
[RMT 1090]. 


P. 224, column 3, formula 1.1.0.1, line 3. For (n — 3) read (nm — 2). 
P. 274, column 3, formula 1.3.1.2, line 2. 


For 
read 


P. 275, column 3, formula 2.3.1.0. For (asin bx + bcos bx) read 
(a cos bx + b sin bx). 


These are in addition to the loose leaf sheet of errata supplied with 
the book. 


MuRLAN S. CORRINGTON 
Radio Corporation of America 
Camden, New Jersey 
224.—NBSMTP, Table of Natural Logarithms, v. 4. New York, 1941. 
Table I, p. 484, in the argument column 


for x = 9.2834 read x = 9.8234 


YVONNE SETTLE 
Midwest Research Inst. 
Kansas City, Missouri 


225.—J. W. NicHotson & J. R. AtrrEy, BAAS Report, 1927. 
In the table of confluent hypergeometric functions, p. 220-254, there 
is an error on p. 243 in M(7/2, 4, .5). 
for x = 1.55306 read 1.55290 


YVONNE SETTLE 
Midwest Research Institute 
Kansas City 2, Missouri 


226.—G. W. SPENCELEY & R. M. SPENCELEY, Smithsonian Elliptic Functions 
Tables, Publ. no. 3863, Washington, 1947. [See MTAC, v. 3, p. 89.] 


The whole E(¢, k) column (p. 71, 73) for the modular angle @ = 18° 
is wrong. It should read: 


r E(¢, k) r E(¢, k) 

0 0.00000 00000 00 

1 0.01789 37526 82 6 0.10732 43169 64 
2 0.03578 64115 08 7 0.12519 52058 38 
3 0.05367 68841 57 8 0.14305 85232 78 
4 0.07156 40813 75 9 0.16091 32180 25 
5 0.08944 69184 96 10 0.17875 82508 48 


is 
4 > 
‘ 
| 


402 E. Church St. 
Oxford, Ohio 
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E(¢, k) 
0.19659 25959 83 
0.21441 52425 30 
0.23222 51958 40 
0.25002 14788 57 
0.26780 31334 34 


0.28556 92216 21 
0.30331 88268 98 
0.32105 10553 91 
0.33876 50370 29 
0.35645 99266 71 


0.37413 49051 83 
0.39178 91804 68 
0.40942 19884 56 
0.42703 25940 37 
0.44462 02919 51 


0.46218 44076 22 


0.47972 42979 52 © 


0.49723 93520 44 
0.51472 89918 93 
0.53219 26730 08 


0.54962 98849 91 
0.56704 01520 56 
0.58442 30334 98 
0.60177 81241 05 
0.61910 50545 24 


0.63640 34915 61 
0.65367 31384 42 
0.67091 37350 12 
0.68812 50578 87 
0.70530 69205 52 


0.72245 91734 12 
0.73958 17037 93 
0.75667 44358 87 
0.77373 73306 65 
0.79077 03857 25 


0.80777 36351 12 
0.82474 71490 80 
0.84169 10338 19 
0.85860 54311 38 
0.87549 05181 05 


E(¢, k) 
0.89234 65066 55 
0.90917 36431 51 
0.92597 22079 17 
0.94274 25147 28 
0.95948 49102 74 


0.97619 97735 85 
0.99288 75154 29 
1.00954 85776 81 
1.02618 34326 61 
1.04279 25824 46 


1.05937 65581 55 
1.07593 59192 16 
1.09247 12526 02 
1.10898 31720 49 
1.12547 23172 55 


1.14193 93530 55 
1.15838 49685 84 
1.17480 98764 16 
1.19121 48116 95 
1.20760 05312 41 


1.22396 78126 53 
1.24031 74533 91 
1.25665 02698 50 
1.27296 70964 18 
1.28926 87845 34 


1.30555 62017 22 
1.32183 02306 31 
1.33809 17680 52 
1.35434 17239 46 
1.37058 10204 46 


1.38681 05908 66 
1.40303 13787 01 
1.41924 43366 20 
1.43545 04254 61 
1.45165 06132 12 


1.46784 58740 01 
1.48403 71870 72 
1.50022 55357 67 
1.51641 19065 02 
1.53259 72877 46 


G. W. SPENCELEY 
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11 51 
12 52 
13 53 
14 54 
15 55 
16 56 
: : 
18 58 
19 59 
20 60 
21 61 
| 22 62 ; 
23 63 
24 64 , 
25 65 
26 66 
27 P| 67 
28 68 
30 70 
31 71 
32 72 
33 73 
34 74 
35 75 
36 76 
37 77 
38 78 
39 79 
40 80 
41 81 
42 82 
43 83 
44 84 
45 85 
46 86 
47 87 
48 88 
49 89 
50 90 
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UNPUBLISHED MATHEMATICAL TABLES 


156[A].—M. K. Pace & M. K. PretL, Further Computation of the Digits of e. 
One page tabulated from punched cards. Deposited in the UMT FILe. 


This gives the value of e to 3333D as obtained on the IBM 604 calculat- 
ing punch. This work supplements that of F. GRUENBERGER and O. Mar- 
LOWE as reported in MTAC, v. 6, p. 123-124. The remainders of the various 
terms at 3333D and 3347D will be furnished upon request. 

M. K. PaGe 

M. K. PFEIL 


Numerical Analysis Laboratory 
University of Wisconsin 
Madison, Wisconsin 


157[A].—J. W. WrReENcH, JR., Log and Related Values. One typewritten 
page deposited in the UMT Fie. 
The values of In z, log x, $ In 2x, $ log 2x, and log 2 are given to 329D. 
The first of these was computed from 


In = In (2"-3-78-175) — 4 log (5%-11-13) + log (1 — 6) 


where ¢ is approximately 2.51869-10-*. These values confirm the 214D 
values published by UHLER.! 


4711 Davenport St., N.W. 
Washington 16, D. C. 


1H. S. UBLER, aeons other basic constants,” Nat. Acad. Sci., Washington, Proc., 
v. 24, 1938, p. 23-30 


158[F].—A. GLoDEN, Solutions of x*+1=0 (mod p) for 600000 < p 
< 800000. Typewritten manuscript, 22 p. Deposited in the UMT FILE. 
This is a continuation of two previous tables for p < 600000 [MTAC, 
v. 3, p. 96]. The present table gives two solutions for every possible prime p 
between 600000 and 800000. 

With the help of these tables various results on the factors of x* + 1 were 
obtained, such as 


J. W. WRENCH, JR. 


8204 + 1 = 626929-721169. 
These are appended. 


11 rue Jean Jaurés A. GLODEN 


Luxembourg 


AUTOMATIC COMPUTING MACHINERY 


Edited by the Staff of the Machine Development Laboratory of the National Bureau 
of Standards. Correspondence regarding the Section should be directed to Dr. E. W. 
CANNON, 415 South Building, National Bureau of Standards, Washington 25, D. C. 


TECHNICAL DEVELOPMENTS 


THE ELECTRONIC COMPUTER AT THE INSTITUTE 
FOR ADVANCED STUDY 


An all-electronic, general purpose, digital computer has been in opera- 
tion at the Institute for Advanced Study in Princeton, New Jersey, since 
June 1952. 
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The overall characteristics of the Institute machine are displayed below 
and may be seen to follow closely the requirements established by Burks, 
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GOLDSTINE and VON NEUMANN in June 1946. 


MACHINE CHARACTERISTICS 


OVERALL ~ 
Binary 
Parallel 
40 bit words 


2300 thermionic tubes 
2’ X 6’ X 8’ in size 
16 kilowatts of power 


MEMORY 


William’s type electrostatic storage 
1024 word capacity 

40 parallel stages 

Standard 5CP1A cathode ray tubes 
Synchronous memory control 

24 u.s. period 


CONTROL 


1-Address code with each 40 digit 
word containing two 20 digit orders 

Asynchronous 

Direct coupled 

20 basic instructions 

10 bit instructions 


ARITHMETIC 


Completely direct coupled 

Fixed binary point 

Complemented representation of 
negative numbers 

All numbers —-1<x<+1 

12 carry time 

60 yu.s. addition order time 

700 u.s. average multiplication order 
time 

1100 yu.s. division order time 


INPUT-OUTPUT 


IBM PUNCHED CARD 


12 words per card 

40 binary digits per word 
100 cards per minute—load 
100 cards per minute—punch 


TELETYPE TAPE 


5 holes across 

4 bits information 

8 minutes per 1024 words—load 

32 minutes per 1024 words—print 
16 minutes per 1024 words—punch 


MAGNETIC DRUM 


In process 


The translation of the preliminary principles of logical design into a 
working high-speed machine was accomplished by a small group of engineers 
and technicians directed by JULIAN H. BIGELow. Bigelow’s design disci- 


plines were: 


(a) Compactness, achieved by seeking optimum geometrical distribution 
of components, to minimize the effect of the purely parisitic inter- 


connections. 


(b) Minimization of the number of tubes, achieved by (a) along with 
the elimination of logical redundancies and special purpose circuitry. 
(c) Logical operations independent of signal waveforms. 
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(d) Sequential control operation with each key operation initiated by 
the safe completion of the preceding operation. 

(e) Safe circuit operation between essentially zero duty cycle and its 
maximum safe repetition rate since the duty cycle over a long period 
of varied computations could not be predicted. 

(f) All of the above bounded by very conservative design criteria for 
the components. 


The major departure of the final machine from the proposals of the 1946 
report is the use of the Williams storage phenomena in a cathode ray tube 
memory. The 40-stage memory with a capacity of 1024 words was developed 
by James H. PoMERENE and Bigelow. 


2 
> 
Digit Resolver 
Adder 
Williams 0 
Tube 0 Controls 
0 Com, ntl 
Gates 
0 y 


Williams 
Memory 


Scnematic of Basic Connections Between Wiutiams 


Memory, Conrror, ann Unirs 


Fic. 1. 
Structure. 


The frontispiece shows a view of the Institute machine. 

Figure 1 is a schematic establishing the physical layout of the machine 
and displaying the interconnections between the Memory, Control, and 
Arithmetic Units. 
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The large block at the bottom contains the 40 cathode ray tubes (SCP1A), 
20 on each side. Each memory tube has an associated three stage amplifier 
and an eight tube discriminator unit which controls the beam turn-on. 

The upper portion contains the arithmetic and control organs. On the 
left there are three double rank registers each containing two rows of 40 
flip-flops and two rows of gate tubes for entry into them. The left-most stage 
contains the most significant digit. The machine is designed to accommodate 
numbers in the range — 1 < x < + 1, with negative numbers handled in 
the complemented representation. The most significant digit thus becomes 
the sign indicator with the binary point fixed to the right of it. 

The register labelled RI is called the Accumulator. It is a double rank 
shifting register which holds the resident number and then the final sum 
during an addition, the partial remainders during division and the most 
significant digits of the partial products during multiplication. 

RII is called the Arithmetic Register. It is a double rank shifting register 
which holds the multiplier and the least significant digits of the partial 
products during multiplication and the quotient during division. 

RIII is a non-shifting Memory Register. In the upper row of 40 flip-flops 
it receives numbers from the memory and presents them to the adder via the 
40 complementing gates. It contains the incident number during addition, 
the multiplicand during multiplication and the divisor during division. In 
the lower row of flip-flops it receives order-pairs from the memory, each half 
containing a ten digit address and a ten digit instruction. Since order-pairs 
and numerical data are indistinguishable outside of the process in which they 
are used, any order-pair may be altered by bringing it from the memory to 
the arithmetic organ via the upper row of RIII, performing an arithmetic 
operation on it and storing the result back in the memory. 

On the right the row of adder chassis contains 40 parallel circuits each 
receiving three inputs; the resident digit from RI, the incident digit from 
RIII (upper) via the complement gates, and a carry digit from the next less 
significant stage. The adder circuit contains four tubes which propagate a 
carry to the next stage and reduce the eight possible inputs to four possible 
voltage levels. The row of chassis above the adder is called the digit resolver 
and reduces the four voltage levels coming out of the adder to two voltage 
levels representing a “1” or a “0”. This information is then communi- 
cated to RI. 

The lower row of chassis comprise the control organ containing four 
major units: the main control, the shift counter, the address counter, and 
the memory control. 


1) The main control establishes the sequence for observing the two 
halves of any order pair and bringing up a new order. It decodes the order 
digits and commands the execution of the resulting process. It is completely 
asynchronous. It does not go through the expense of matrix decoding, but 
rather utilizes unique logical combinations of several order digits to com- 
mand the desired process. 

2) The shift counter counts the number of shifts in the shifting registers 
and contains recognition circuits which determine when shifting shall be 
terminated appropriately to the order being done. 


by 

its 

od 

‘or 

46 

be ‘3 

ed 

ie 


112 AUTOMATIC COMPUTING MACHINERY 


3) The memory control with its associated pulse chain establishes read- 
ing and writing routines under the command of the main control, and, 
independently, routines to control the regeneration of the static information 
in the memory. 

4) The address counter remembers the locations of the last point re- 
generated on the face of the cathode ray tube and of the last order-pair 
executed; advancing these addresses by one each time. It is also capable of 
receiving from RIII, however, the address associated with any instruction. 
During every memory cycle one of the above addresses is used to command 
the deflection of the cathode ray beams. 


Memory. 


During any period of time in which the memory is not specifically in- 
volved in delivering information to or receiving information from the 
arithmetic organ, the memory control reverts back to sequential regenera- 
tion of the 32 X 32 words stored on the phosphor surfaces of the 40 cathode 
ray tubes. The timing of this process is governed by a clock and associated 
pulse chain, with a 24 u.s. period. Within this period there is one moment 
when the memory control is receptive to intervention by the main control. 
The main control may at this time be specifying one of the following 
processes: 


(a) Bring the next order-pair of a consecutive series of orders into the 
lower row of 40 toggles in the memory register RIII. 

(b) Depart from this consecutive series of orders to bring a new order- 
pair into the lower row of the memory register RIII from an address 
specified by a ‘‘transfer control” order in RIII. 

(c) Bring a number into the upper row of the memory register RIII from 
the address specified by the order in RIII. 

(d) Store the number residing in the accumulator register RI at the 
address specified by the order in RIII. 


The memory control reports back to the main control when it accepts 
the main control intervention and when it completes any of the actions 
described in (a)—(d). In this manner the synchronous memory operation is 
safely incorporated into the asynchronous arithmetic process. 

The inherent coupling between adjacent storage locations on the face 
of any cathode ray tube places a limit on the number of times any point 
on the raster may be consulted before its neighbors are regenerated. For 
the Institute machine this number, called the ‘‘read around” ratio, is in the 
neighborhood of 30. 


Arithmetic and Control Operations. 


The twenty basic orders executed by the Institute machine are listed 
below with the average time per order measured in a long sequence of 
similar orders. Many of the processes are operating with inordinate factors 
of safety. Thus the speed figures below reflect the present arbitrary state of 
adjustment of the machine, not its maximum capabilities. 
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Bring a number from position x in the memory and add it to the number 
residing in the accumulator with the following choices: 


(a) “Clear” or “Hold” where “‘Clear’’ signifies the desire to preclear 
the accumulator to zero. 

(b) “Magnitude” or “Number” where ‘“Magnitude’”’ signifies that the 
absolute magnitude of the quantity at x is to be added. 


(9) Multiply: Multiplier all ones... 990 y.s. 


Bring a number from position x in the memory and multiply it by the 
number residing in the Arithmetic Shifting Register, RII. This order has 
choice of ‘Round off.” 


Divide the number in the Accumulator by a number to be brought up 
from position x in the memory. 


Bring a number from position x in the memory and place that number 
in the Arithmetic Shifting Register, RII. 


42 us. + 8 w.s./shift 

Shift the numbers in the Accumulator register, RI and the Arithmetic 
Shifting Register RII, »(< 47) places to the left. During this process the 
contents of RI shift end around into the right end of RII. However, only 
zeros are propagated through the right end of RI. The number 1 is stated 
in the address portion of the order. 


Shift the numbers in the Accumulator register RI and the Arithmetic 
Shifting Register RII n(< 47) places to the right. During this process the 
contents of RI shift end around into the left end of RII. However, only the 
sign digit of the original quantity in RI is propagated through the left end 
of RI. The right shift order has the choice of Round Off. If Round Off is 
chosen, 2—** is added into the quantity in RI at the end of ‘‘n”’ right shifts 
and a total of m + 1 shifts are executed. 


Take the number existing in the Accumulator register RI and store it 
at position x in the memory. 


Bring a new order pair into the lower rank of the Memory Register, 
RIII from position x in the memory and proceed to observe that order which 
occupies the same half of the order pair as the unconditional transfer being 
executed. 


d- 

on : 

e- 5 

ir 

of 

n. 

id 

le 

d 

it 

1. 

ig 

n 

e 

s 

Ss 

e 

t 

r 


114 AUTOMATIC COMPUTING MACHINERY 


(16) Non-Parity Unconditional Transfer... 39 us. 
Bring a new order pair into the lower rank of the Memory Register, 
RIII from position x in the memory and proceed to observe that order which 
occupies the half of the order pair opposite to that occupied by the uncondi- 
tional transfer being executed. 


(17-18) Parity and Non-Parity Conditional Transfers.............. 39 yu.s. 
If the number in the accumulator is positive, execute operations de- 
scribed in 15 and 16. 

If the number in the accumulator is negative, ignore the request for 
transfer and proceed to the next order in sequence. 


Take the number in the Arithmetic Shifting Register RII and add it to 
the quantity in the accumulator register RI. This order includes all of the 
eight variants described in orders (1-8). 
(20) Input-Output Order 
(a) IBM 
Transfer 12 » words from » IBM punched cards to 
the 12 » memory addresses starting at address x. 
Transfer 12 » words from the memory addresses 
starting at address x to m IBM cards. 
(b) Magnetic Drum 


Transfer 32 words starting with block mp of Track A 
or B from the magnetic drum to the 32 m memory ad- 
dresses starting at address x (1<n<32; 0<mo<31). 
Approx (150 + 3) milliseconds. 
Transfer 32” words from the memory addresses 
starting at address x to the m Blocks of track A or 
B of the magnetic drum following the starting 
Block no. 


Operating Experiences. 

During the first year of machine operation, a considerable amount of 
time was spent in strengthening weak points in the machine. Towards the 
latter part of the year the machine came under control with a routine main- 
tenance program permitting long periods of uninterrupted computation. 
The bulk of the machine’s useful computation time has been devoted to 
the solution of sets of partial differential equations representing simplified 
models of the atmosphere. These problems were formulated by a group led 
by J. G. CHARNEY. In addition, problems in number theory, matrix inversion, 
and integral equations have been solved by its circuitry. These problems 
required about 3 X 10*® multiplications with the time split fairly evenly 
between automatic computation and input-output. 

GERALD EsTRIN 
Institute for Advanced Study 

Princeton, N. J. 
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1016. BRUNO DE FinettT1, Macchine “che pensano” (e che fanno pensare), 
Pubblicazioni delle Facolta di Scienze e di Ingegneria dell’Universita di 
Trieste, 1952, 34 p., 23.4 X 29.5 cm. 


This brochure contains an expository treatment of the automatic digital 
computer of our era, seen through the eyes of a mathematician and philoso- 
pher. Technical details of machine design, being of minor importance to 
him, are relegated to an appendix which includes photographs and descrip- 
tions of some individual machines and the inevitable summary table of 
speeds and memory capacities. While this appendix is purposely incomplete 
and not too, stimulating, the main part of the booklet is written with a depth 
of understanding and a clarity and conciseness of expression seldom found 
in this field. To a reader who knows nothing of computing machines, it 
explains the relationships between their development and such fields as 
neurophysiology and psychology, formal logic, information theory, and 
cybernetics. 

It explains the basic functions and components of automatic digital 
computers, the essential steps involved in programming and coding, the 
design of “subroutines,” and their assembly into comprehensive routines. 
It touches upon the major classes of computational problems, evaluation of 
functions and use of tables, matrix inversion and numerical integration of 
ordinary and partial differential equations. There is a section on relaxation 
methods and one on the Monte Carlo method. In each of these subjects the 
fundamental features are adequately brought out in just a few sentences. 
While it is not clear why the author arranged his topics the way he did, their 
selection and presentation leaves nothing to be desired. 


F. L. Att 
NBSNAML 


1017. F. GRUENBERGER, Computing Manual, The University of Wisconsin 
Press, Madison, 1952, 123 p., 15.5 X 23.5 cm. Price $2.00. 


This manual was used as a text in a course entitled ‘Theory and Opera- 
tion of Computing Machines” held at the University of Wisconsin during 
1951 and 1952. The author discusses IBM techniques for chi-squared 
analysis, random sampling, questionnaires, bivariate tables and differercing. 
He also presents miscellaneous operating hints for the IBM laboratory 
worker. 


R. K. ANDERSON 
NBSCL 


1018. H. D. & V. R. Huskey, “Electronic computers aiding management 
control,” The Journal of Accountancy, Jan. 1952, p. 69-75. 


In this very informative article, the authors point out the many econo- 
mies which may be effected by the use of electronic computation in business 
administration. Not only would present-day calculations be carried out 
more efficiently, but the savings in labor, time, and expense would enable 
various organizations to undertake the solutions of problems hitherto con- 
sidered prohibitive, for example, the forecasting of trends, as based upon ac- 
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cumulated data. The recent election prediction carried out by the UNIVAC 
amply confirms this statement made several months before Election Day 
by the authors. 

The authors indicate some of the factors which account for the economies 
realized from the use of the electronic computers, such as the immense speed 
of internal computation. For example, in the time it takes to say the word 
“multiplication,” some computers can perform about 2500 multiplications, 
or 6000 additions, of pairs of ten-digit numbers. Moreover, these machines 
are constructed to recognize, and cause automatically the execution of, 
instructions to perform numerical operations. The resulting elimination of 
human intervention is a potent time-saver indeed; it also permits a consider- 
able reduction in personnel and in space requirements. 

A brief but illuminating history of high-speed computing machines con- 
cludes this excellent article. 

The reviewer is left with the feeling, however, that the authors give too 
optimistic a concept of the prowess of the present-day machines as applied 
to business needs. The prodigious feats of computation which they have 
exhibited do not apply to all classes of problems. Where defense tasks are 
concerned, the astounding reduction in time required for their solution is 
an advantage far outweighing any other considerations. But business 
managers cannot expect any spectacular savings in dollars and cents until 
efficient electronic sorter-collators and files allowing random access have 
been perfected. 


IpA RHODES 
NBSCL 


1019. H. KENosIAN, “Unitized pulse circuits speed computer design,” 
Electronics, v. 25, no. 10, Oct. 1952, p. 156-157. 


These unitized pulse circuits were originally developed in Project 
Whirlwind at MIT. They have been further developed and made commer- 
cially available by the Burroughs Adding Machine Co. The ten types of 
unitized circuits include generators for half-sine-wave pulses of 0.1 micro- 
second duration, a flip-flop which can be triggered by the pulses and gives 
out d-c levels for gating, detectors of coincidence between a pulse and a d-c 
level or between up to five sources of d-c levels, pulse delays, pulse mixers, 
and channel selectors. Each unit occupies 3} inches of height in a 19 inch 
relay rack. They may be interconnected by 93 ohm coaxial cable into com- 
plex systems. For example, the Burroughs Laboratory Computer incor- 
porates nearly 500 of these basic pulse control units. 


R. D. ELBouRN 
NBSEC 


1020. H. Koppet, “Digital computer plays NIM,” Electronics, v. 25, Nov. 
1952, p. 155-157. 


Four stacks of seven chips each are indicated by lights. A “‘safe’’ starting 
condition makes it possible for the player to beat the machine, but not by 
simply memorizing plays. No electrical diagram appears, but a block dia- 
gram shows adders, binary counters, gates, pulser, and scanner. Other 
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binary game machines at University of Toronto, NBS Corona Labs, etc., 
are noted. The win theory of this game is given by BouTon.'! 

J. Howarp WriGcHT 
NBSEC 


1C. L. Bouton, “Nim, a game with a complete mathematical theory,” Annals of Math., 
s. 2, v. 3, p. 35-39. 


1021. F. C. MuLianey, “Design features of the Era 1101 Computer,” 
Electrical Engineering, Nov. 1952, p. 1015-1018. 


The article gives a brief description of a computer using a magnetic 
drum memory. The first part deals with the basic facts of the machine: what 
it consists of, what it does, and how it does it. The second part is devoted 
to the constructional features of the machine and is provided with illustra- 
tions. The third section treats the problem of testing and maintenance. 
Finally, the operational record for over 7000 hrs. is given and a resume of 
desirable features of this computer presented. 

A. WEINBERGER 
NBSEC 


1022. OrFicE oF NAVAL RESEARCH, Digital Computer Newsletter, v. 4, no. 4, 
Oct. 1952, 11 pages. 


The present status of the following digital computers is treated briefly 
in this number: 


1. Naval Proving Ground Calculators 
2. The UNIVAC 
3. Whirlwind I 
4. Computer Research Corporation Computers 
CADAC 102-A 
CRC 105 
CRC 107 
Moore School Automatic Computer (MSAC) 
. Raytheon Automatic Computer (RAYDAC) 
. The SWAC 
. Aberdeen Proving Ground Computers 
The ORDVAC 
The EDVAC 
The ENIAC 
The Bell 
IBM 
9. C.S.I.R.0O. Mark I 


Under Data Processing and Conversion Equipment: 
1. Graph Plotter 


1023. L. PacKER & W. J. Wray, JR., “Germanium photodiodes read com- 
puter tapes,” Electronics, v. 25, Nov. 1952, p. 150-151. 


A six-channel device for reading paper tape similar to the teletype device 
uses germanium photoelectric sensing at a speed of 100 pulses per second, 
although the electronic circuits are capable of 35,000 per second. Small 
size and large output of the germanium units are advantageous. Only six 
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diodes are used, but nonetheless this reviewer suggests that reliability is 
doubtful on the basis of various reports of instabilities due to the plastic 
encapsulation of these particular diodes. 


J. HowarD WRIGHT 
NBSEC 


1024. J. RaBinow, ‘‘The notched-disk memory,” Electrical Engineering, v. 
71, 1952, p. 745-749. 


One of the increasingly pressing problems of government and business 
is that of the compact storage of large amounts of information in a quickly 
accessible form. 

This paper describes a storage device consisting of many thin iron oxide 
coated aluminum disks. Information may be recorded on both sides of the 
disks as coded magnetic pulses. The disks are arranged in the form of a 
toroid. When it is desired to read or record information on a particular disk, 
the magnetic recording heads are automatically swung to the chosen disk 
and the disk is rotated between the recording heads for one revolution. 

It is estimated that such a device consisting of 588 disks 20 inches in 
diameter would have a capacity of a quarter of a billion binary digits. Test 
results demonstrate that stored data can be reached in times of the order 
of half a second. 


E. F. AINSWORTH 
NBSEC 


NEws 


Joint Computer Conference.—On December 10-12, 1952, a second joint annual com- 
puter conference and exhibition was held in New York under the auspices of the American 
Institute of Electrical Engineers, Institute of Radio Engineers, and Association for Com- 
puting Machinery. The program for the meeting was as follows: 


Wednesday, December 10, 1952 Registration and Exhibits 
Morning session S. B. Wittiams, Chairman, Burroughs Add- 
ing Machine Co., Consultant 
Keynote address N. H. Taytor, MIT 


Recording techniques for digital coded A.W. TyLer, Eastman Kodak Co. 
data 
Afternoon session M. M. Astranan, IBM, Chairman 
Converters for teletype tapetoIBMcards G.F. IBM 
Punched card to magnetic tape converter E. BLUMENTHAL and F. Lopey, Eckert- 
for UNIVAC Mauchly Div., Remington Rand, Inc. 
Devices for transporting the recording R.L. Snyper, Jr., Consulting Engineer 
media 
Buffering between input-output and the A. L. Lerner, NBS 
computer 


Thursday, December 11, 1952 


Morning session W. S. MacWit1aMs, Chairman, BTL 
SEAC input-output system: 
SEAC input-output system S. GRuNWALD, NBS 


Input-output devices used with SEAC J. L. Pike 

Auxiliary equipment to SEAC input- Ruta C. HavetTer, NBS 
output 

Operating experience E. AinswortH, NBS 
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ity is UNIVAC input-output system: 
lastic The Uniservo-tape reader and recorder H. F. WetsH and H. Luxorr, Eckert- 
Mauchly Div., Remington Rand, Inc. 
IT Input devices L. D. Witson and E. RoGGENSTEIN, Eckert- 
Mauchly Div., Remington Rand, Inc. 
Output devices E. Masterson and L. D. Witson, Eckert- 
ng, Vv Mauchly Div., Remington Rand, Inc. 
' Afternoon session J. H. Howarp, Chairman, Burroughs Add- 
ing Machine Co. 
siness RAYDAC input-output system: 
lickly RAYDAC input-output systems W. Gray and K. Ragier, Raytheon 
System design L. Fern, Raytheon 
oxide Operating experience with RAYDAC F. Dean, Raytheon 
f the IBM 701 input-output system: 


Type 701 input-output organization L. D. Stevens, IBM 


of a Type 726 magnetic tape reader and W. S. Busuix, IBM 
disk, recorder 
disk Magnetic tape techniques and per- H.W. NorpyKe, IBM 
formance 
es in Friday, December 12, 1952 
Test Morning session J. C. McPuerson, Chairman, IBM 
order Survey of analog-to-digital data con- H. E. Burke, Consolidated Engineering 
verters Corp. 
7 Survey of mechanical type printers J. Hosxen, A. D. Little, Inc. 
Survey of nonmechanical type printers Lt. R. J. Rossnem, GWU 
Anelex printer L. Rosen, Anderson-Nichols Co. 
Eastman printer R. G. Taompson, Eastman Kodak Co. 
Afternoon session J. G. Bratnerp, Chairman, Moore School of 
onal Electrical Engineering, Univ. of Penn. 
Pore. Ferranti input-output equipment B. W. Po.varp, Ferranti Ltd., London 
_ Garment tag equipment O. G. Hesse, Sears Roebuck and Co. 
Numerically controlled machine tool W. Pease, MIT 
Summary and forecast S. N. ALEXANDER, NBS 
Add- Informal discussion 
OTHER AIDS TO COMPUTATION 
BIBLIOGRAPHY Z 
1025. LEonrLpA ALTMAN, “An analog computer as an aid to shower theory 
“2 calculations,” Rev. Sci. Instruments, v. 23, 1952, p. 382. 
. A. c. resolvers are used to solve certain trigonometrical equations. 


F. J. M. 


1026. S. Boswortu, “A new analog computer,” Electronics, v. 24, no. 8, 
Aug. 1951, p. 216-224. 


This article describes ‘‘an inexpensive desk-size electronic differential 
analyzer’’ suitable for differential equations with constant coefficients. The 
device contains 20 d.c. amplifiers, 23 ten-turn potentiometers, and 8 inte- 
grating capacitors. The problem is set up on terminal boards which are 
plugged into the device. 


F. J. M. 


. 
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1027. F. E. Brooks, Jr. & H. W. Smitu, “A computer for correlation func- 
tions,” Rev. Sci. Instruments, v. 23, 1952, p. 121-126. 


A continuous computer is described for calculating fo'(F(#)G(t + 7) 


— FG(t + 7))dt for a given value of 7, where F(t) and G(t) are given arbi- 
trary functions of ¢ and F is the average value of F on the given interval of 
integration. The correlation function is computed from a number of such 
expressions. F(t) and G(¢) are recorded on magnetic tapes in the form of pulse 
length modulation of a rectangular wave with a repetition rate of five cycles 
per second. The mathematical operation required is performed by a multi- 
plier integrator. The integration is performed by a condenser which is fed 
by a multiplying unit whose output comes through a gate, which is controlled 
to be off or on by the pulses for one factor. The input to this gate is the result 
of demodulating the signal for the other factor and obtaining an amplitude 
modulated signal. The variable r is introduced by varying the relative posi- 
tion of the two tapes. It is planned to use this in the study of the fading of 
radio signals due to various factors and other atmospheric studies. An 
example involving thermistors is given. 
F. J. M. 


1028. F. B. Danrets, “A planimetric method of harmonic analysis,”’ Rev. 
Sci. Instruments, v. 23, 1952, p. 369-370. 


This article gives a formula in terms of the coefficients of the Fourier 
series of f(x) for an integral which is readily evaluated by a planimeter. 
Let gx(x) = (— if j/2ke < x < (j + 1)/2km. The integral obtained by 
the planimeter is /<?*¢,,fdx. This can be considered as an approximation to 
a, the coefficient of cos kx in the Fourier series for f(x). A phase shift yields 
an approximation to };, the coefficient of sin kx. 

F. J. M. 


1029. E. A. GoLpBERG, “Step multiplier in guided missile computer,” 
Electronics, v. 24, no. 8, August 1951, p. 120-124. 


A set of relays can be used to obtain a conductance whose value is 
determined by a binary coded number .aiaz2 --- a». Thus if a; = 1, the 
jth relay contact will introduce a conductance with value 2’ units. Now 
suppose the state of the j’th relay is determined by the (m — j)’th stage of 
a binary counter. Apply a fixed voltage to this conductance and compare 
the resulting current with another current whose value corresponds to a 
mathematical variable x. The difference controls gates across the input of 
the binary counter, so that if there is a significant discrepancy between these 
currents, an appropriate gate is opened and counting pulses reach the 
counter causing it to count up or down until the value of the conductance 
corresponds to the value of the variable x. The same binary counter can be 
used to control any number of such relay conductances and thus any number 
of conductances with the value of a specified mathematical variable x can 
be produced. If now, voltages corresponding to other variables, y:, ---, yr, 
are applied to these conductances, currents with values xy, ---, xy, may 
be obtained. 

The actual conductances are. T networks and 11 stages are used. The 
counting pulses are at a rate of 1 kilocycle and the relay conductances in- 
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volve a specially developed relay with a response time of approximately 100 
microseconds. The circuit for the binary counter is also discussed. 

These multipliers were developed as part of “‘Project Typhoon” for the 
U. S. Navy. They are used in the Typhoon guided missile simulator where 
multipliers are needed with a precision higher than that available from servo 
multipliers. The article also discusses various features of this computer, 
which is essentially a large electronic differential analyzer using 445 stabilized 
d.c. amplifiers. These can be used for either integration or summation. There 
are 19 computing servo units and 36 multiplications can be made by the 
above described stepmultipliers (involving the 9 distinct x factors). 

F. J. M. 


1030. A. E. Hastincs & J. E. Meape, “‘A device for computing correlation 
functions,” Rev. Sci. Instruments, v. 23, 1952, p. 347-349. 


This device produces the integral /_r7 f(#) f(t + r)dt. The function f(é) 
is recorded on a magnetic tape by amplitude modulation and read by two 
heads whose separation corresponds to r. One head signal is used to produce 
a current, the other a voltage and the result is obtained from a wattmeter. 

F. J. M. 


1031. INSTITUTE oF Gas TECHNOLOGY, ‘‘Gas distribution system calculator,” 
Rev. Sci. Instruments, v. 23, 1952, p. 61. 


A photograph of a special purpose network analog computer and a short 
description are given. 
F. J. M. 


1032. R. D. Minpiin & M. G. Satvaport, “‘Analogies,’”’ Handbook of Ex- 
perimental Stress Analysis. John Wiley & Son, New York, 1950, Chapter 
16, p. 700-827. 


When the behavior of two physical systems can be described in identical 
mathematical form, they are said to be analogous. If experiments and meas- 
urements can be made with greater ease and accuracy on one of the systems, 
then it may be used as a computing device for the solution of complex prob- 
lems arising in the more intractable system. In this article, the authors 
present a comprehensive and detailed survey of those analogies between 
elasticity and other physical phenomena which have been used in the quan- 
titative determination of states of stress in elastic bodies. Each analogy is 
treated carefully, deriving first the mathematical equations governing both 
the original physical system and the analogous system. The corresponding 
quantities in the two systems and the required scale factors are then readily 
recognized. Methods of physically realizing and measuring quantities in the 
analogous system are then described, including a discussion of the overall 
accuracy obtainable. An excellent bibliography containing sixty-seven 
entries is appended. 

The earliest use of analogies in elasticity theory was apparently in the 
determination of stresses in a cylindrical bar bent by a terminal load or 
twisted by terminal couples. At a sufficient distance from the end of the 
bar, the problem may be reduced to the solution of Laplace’s or Poisson’s 
equation with suitable boundary conditions. Therefore, any of the other 
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physical phenomena satisfying these equations may be used as analogies. 
PRANDTL, in 1903, suggested the use of a homogeneous membrane, under 
uniform tension per unit of length and subjected to a normal pressure per 
unit of surface area, the deflection of which satisfies Poisson’s equation, 
provided the slope of the membrane relative to a plane is everywhere small. 
This analogy has been extensively employed and a full description is given 
of methods of its physical realization, of fixing its boundaries so that the 
boundary conditions are analogous to those of the original problem and 
methods of measuring deflections of the membrane accurately. Soap films 
have been used primarily, although the use of rubber and the separation 
surface of two immiscible liquids are described. The steady-state electrical 
potential in a thin plate of constant thickness, which satisfies Laplace’s 
equation, has also been utilized. 

The use in analogies of continuous physical systems such as the mem- 
brane or thin conducting plate, as opposed to lumped parameter systems, 
has the advantage that the analogous system satisfies the same system of 
equations as the original system, whereas the lumped parameter system will 
satisfy only a finite difference approximation to those equations. However, 
the use of continuous systems requires usually the construction of new 
apparatus (molds, etc.) for the solution of each problem and delicate experi- 
mental technique for accurate measurements. For these reasons, it appears 
to the reviewer, that a lumped D. C. network would represent a more flexible 
analogy to the torsion problem than the membranes or conducting plates 
heretofore used. 

Analogies for the solution of many other problems in elasticity are 
discussed. Among them are elgctrical analogies for the torsion of bars of 
varying circular section and analogies between problems in elasticity which 
are difficult to study experimentally and other elasticity problems which are 
more susceptible to direct measurement. Thus an analogy exists between the 
stresses produced by a two dimensional steady-state temperature distribu- 
tion and the stresses produced by dislocations, and the latter problem is 
easier to treat experimentally. 

Considerable space is devoted to a lucid exposition of the A. C. electric 
network analogue devised by G. Kron for the general case of static plane 
stress. A lumped parameter system, utilizing inductances and capacitances, 
it satisfies a finite difference approximation to the equations governing the 
continuous system. Its great generality is very attractive. However, the 
analogy requires that the inductances and capacitances have no resistive 
losses. In an actual realization of the analogue described in this article, 
resistive losses resulted in large errors. 

J. H. WEINER 


Columbia University 
New York 27, New York 


1033. REEVES INSTRUMENT CORPORATION. Project Cyclone Symposium II 
on Simulation and Computing Techniques. New York City, April 28- 
May 2, 1952, 237 p. 


The conference reviewed below was the second symposium on REAC 
techniques sponsored by the Reeves Instrument Corp. The first symposium 
was held in 1951 [MTAC, v. 5, p. 250-253]. Project Cycione at Reeves is 
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sponsored by the U. S. Navy Bureau of Aeronautics and the U. S. Navy 
Special Devices Center. The functions of Project Cyclone are: “(1) the 
development and test operation of a guided missile simulator . . . ; (2) the 
operation of a computer laboratory where classified problems in other fields 
may be studied and analyzed; (3) the investigation of new applications for 
Cyclone equipment; (4) theoretical research to ascertain and improve the 
reliability and accuracy of computer solutions; and (5) the coordination of 
the experiences of individual computer groups . és 

The 24 papers presented in this publication can n be roughly classified into 
four types. (a) Papers dealing with circuitry, techniques and auxiliary equip- 
ment to improve the performance of REAC equipment; (b) Papers dealing 
with the application of REAC equipment to various physical problems en- 
countered by the various users of the equipment; (c) Papers dealing with 
computers other than REAC; (d) Papers dealing with basic theory and 
new equipment. 


(a) Circuitry 


Paper No. 10, “Some REAC techniques employed at the David Taylor 
Model Basin’’ by L. Pope discusses techniques for generating input func- 
tions and the use of the REAC as.an harmonic analyzer in problems that 
arose at the DTMB. Paper No. 11, ‘Simulation studies of a relay servo- 
mechanism” by N. P. ToMLINSON considers methods for the analog study 
of a non-linear relay (that is, sampling) servomechanism. Paper No. 13, 
“Precision in high-speed electronic differential analyzers’’ by H. BELL, Jr. 
& V. C. RIDEOUT is concerned with solving linear differential equations with 
variable coefficients and non-linear equations by analog computers and a 
comparison of the accuracy of solutions compared with digital solution. 
Paper No. 15, ‘‘Solution of linear differential equations with time-varying 
coefficients by the electronic differential analyzer’”’ by C. E. Howe, R. M. 
Howe & L. L. Raucs discusses techniques for solving linear differential 
equations on the REAC. Examples include Bessel’s equation, Legendre’s 
equation and the deflection equation [EI(x)y’’] — u(x)d\*y = 0 of simple 
beam theory. Paper No. 19 is entitled “Automatic REAC operation for 
statistical studies” by R. R. BENNETT & A. S. FuLton. Determination of 
the response of a non-linear system to statistical inputs can most conveni- 
ently be carried out experimentally, that is, on an analog computer. Methods 
for changing a Gaussian distribution by passing the noise through a function 
generator are discussed. Paper No. 21, ‘(Checking analogue computer solu- 
tions’”’ by W. F. RicHmMonp, Jr. & B. D. LOVEMAN considers procedures to 
verify analog computer solutions of linear and non-linear differential equa- 
tions. Linear equations are treated by frequency (steady state response) or 
selected root method (transient response). For non-linear systems “dynamic 
substitution”’ is used. It consists in breaking-up the time interval into small 
pieces and substituting the solution on each small sub-interval into the 
original differential equation. Paper No. 24, ‘Modifications and additions 
to the REAC”’ by J. W. Fo tin, Jr., G. F. Emcn & F. M. WALTERs dis- 
cusses circuitry details and auxiliary equipment for improving the operation 
of the REAC. These include a low frequency noise generator, a programmer 
for turning the REAC and recorder on and off, air conditioning the racks, 
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power supply changes, circuit changes in amplifiers and techniques of inter- 
connection of machines. Paper No. 25, ‘‘REAC servo response” by A. H. 
MILLER discusses methods for introducing variable coefficients and non- 
linearities into REAC equipment by the use of servos. A study of the dy- 
namic characteristics of servos for use in the above connection is made. 
Paper No. 26 is ‘“‘Applications of differential relays to solution of REAC 
problems” by L. M. WarsHawsky & W. Braun. A differential relay is a 
single-pole double-throw relay at the output of an operational amplifier. 
Its uses include providing a time-controlled series of switching operations 
and a simulation of ‘‘dead” zone where the operating limits are variable. 
Paper No. 27, ‘“Techniques applied to the design and test of linear filters’ 
by G. NEsTor is concerned with problems met in building and testing a 
Wiener filter. Paper No. 28, ‘‘The role of diodes in an electronic differential 
analyzer” by C. D. Morritt & R. V. Baum discusses the use of diodes to 
simulate certain types of non-linearities. 


(b) Applications 


Paper No. 7: ‘“REAC solution of problems in structural dynamics” by 
C. W. BRENNER. In aircraft design one can no longer consider the air-frame 
as a rigid body. Many components designed only for maximum static loads 
fail in service due to transient loading. This paper considers a simplified 
associated problem of a simply-supported beam struck by a pressure wave 
(a) Normal to the beam, (b) Parallel to the beam. For case (a) there arises 
the differential equation + k’r4¢, = — t/t), r = 1, 2, ---, where 
¢, is the normalized displacement of the rth mode of vibration and k is a 
parameter. A similar equation holds for case (b). These equations are solved 
on the REAC. A discussion of the circuitry details and the troubles that 
were encountered together with their solutions, is presented. 

Paper No. 8: “‘The use of an analog computer and feedback theory for 
the solution of structural problems in the static case” by G. J. MARTIN & 
L. M. Lecatsxr. Structural frameworks (pinned frames, statically inde- 
terminate structures) give rise to certain systems of linear algebraic equa- 
tions of the form 5-71 a:; = M;, i = 1, ---,m, where the 6’s are angles 
representing the rotations of various members of the structures and the M’s 
are moments. These equations bear great similarity to usual mesh equations 
in electrical network theory and are solved on the REAC using a network 
analog. 

Paper No. 9: “Application of the electronic differential analyzer to 
eigenvalue problems” by G. M. Corcos, R. M. Howe, L. L. Rauch & 
J. R. SELLARs. Determination of eigenvalues for Sturm-Liouville systems, 
fourth order equations and partial differential equations. The partial differ- 
ential equations are reduced to ordinary differential equations by separation 
of variables. Since, in general, there are two-point boundary conditions, the 
method is the familiar technique of satisfying the initial conditions at one 
point and varying the solution until the boundary conditions at the second 
point are satisfied. The authors’ experience leads them to conclude that one 
or two end conditions are easy to handle, three end conditions give consider- 
able trouble and four or more end conditions are apt to be hopeless. Four 
examples are considered. 
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Paper No. 12: “Use of the REAC as a curve fitting device” by C. H. 
Murpay. In certain problems in aerodynamics it is necessary to determine 
the coefficients of a differential equation from discrete points on its solution 
curve. The use of the REAC in this connection is described, based on certain 
linearized differential equations of yawing motion. 

Paper No. 14: “An application of analogue computers to problems of sta- 
tistical analysis” by H. H. LANING, Jr. and R. H. Battin. The response x(#) 
of a linear system toa forcing function f(#) isgiven by x(t) = W(t, r) f(r)dr 
where W(t, 7) is the impulsive response (Green’s function, weighting func- 
tion) of the system. A method for computing this quantity is discussed using 
the concept of the adjoint system. 

Paper No. 23: “Solution of partial integral-differential equations of elec- 
tron dynamics using analogue computers with storage devices” by C. C. 
Wane. In certain problems of electron dynamics there arise simultaneous 
partial integral differential equations involving two independent variables. 
How to adapt such problems to REAC analog computers is discussed. 


(c) Non-REAC computers 


Paper No. 16: “JAINCOMP computers and their application to simula- 
tion problems” by D. H. Jacoss. Description of the uses of the JAcobs- 
INstrumentCOMPany digital computer. It is a compact (17 X 21 X 30 
inches), light (110 Ibs.) electronic digital computer. The operations of the 
machine are very rapid (adds two 40 binary digits in eight microsec). Paper 
No. 17: “The decimal digital differential analyzer CRC-105 as a tool for 
simulation and checking analogue computer solutions” by E. Wetss. De- 
scription of the DDDA and its advantages and disadvantages as compared 
with electronic analogue computers and mechanical computers. It operates 
on a decimal number system, has 60 integrators and a maximum accuracy of 
six digits plus sign in each integrator. There are twelve input and twelve 
output channels. 

Paper No. 18: “Problems encountered in the operation of the M. I. T. 
flight simulator’ by W. W. SEIFERT & H. Jacoss, JR. The title accurately 
describes the paper. The simulator is located at the DACL (Dynamic 
Analysis and Control Laboratory) at M. I. T. Paper No. 22: “Analogue com- 
putation of blade designs’’ by D. B. BREEDoN, M. M. Mattuews & E. L. 
HARDER. Calculation of natural frequencies and stresses in vibrating beams 
by the Westinghouse analog computer (Anacom) and checks made on a digi- 
tal machine (IBM Calculating Punch, Type 602-A). 


(d) Theory 


Paper No. 20: ‘Mathematical error analysis for continuous computers” 
by F. J. Murray. This paper presents certain aspects of the work done by 
the author and the reviewer on the mathematical foundations of an error 
theory applicable to mathematical machines. The complete results will ap- 
pear elsewhere. 

Paper No. 29: “A high accuracy time division multiplier” by E. A. 
Go.LpBERG. If a train of rectangular pulses is generated, then the DC com- 
ponent is proportional to the amplitude of the pulses multiplied by the duty 
cycle. This principle is not new. However, by ingenious circuitry including 
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the use of a feed-back system for establishing accurate timing and using a 
current switch essentially independent of tube characteristics, an accurate 
multiplier is obtained which should be useful in electronic analog equipment. 
Using precision resistors, stable voltage sources and steep switching pulses 
an accuracy of about + 0.01% of full scale (75 volts) for a maximum prf 
of 2000 cps is obtainable. 

Paper No. 30: ‘An AM-FM electronic analog multiplier” by W. A. 
McCoo. A signal having both amplitude modulation and frequency modu- 
lation is applied to an (AM) detector and a (FM) discriminator. If it is 
desired to multiply the two quantities X and Y, then the output of the 
discriminator can be made to be X = Y/,Af, and if the output of the detector 
is forced to match the X input plus the reference voltage Yo, then its output 
isZ’ = (Yo + Y)AfandZ = Z’ — X = Yo "XY. (The frequency deviation 
Af is the same for both detectors.) This principle has been used before. By 
the use of feed-back high gain amplifiers, tube characteristics have only a 
secondary effect. A mathematical analysis is made of the circuit showing 
the dependence of Z on the circuit parameters. An experimental model is 
in the process of being developed at the Naval Research Laboratory. 

K. S. MILLER 
New York University 
New York, New York 


1034. DENMAN SHAW & R. N. WELTMANN, “A logarithmic counting-rate 
recorder,”’ Rev. Sci. Instruments, v. 23, 1952, p. 528-531. 


This device uses a spiral cut on a cylinder to record the logarithm of the 
time for a specified number of counts from a Geiger counter and scaler. 
F. J. M. 


1035. S. Sorrky & J. JUNGERMAN, “Electrolytic tank measurements in 
three dimensions,’’ Rev. Sci. Instruments, v. 23, 1952, p. 306-307. 


The authors indicate the techniques necessary to obtain good results 
in the three dimensional case and obtain in a specific example, an accuracy 
of .2% in the potential and 2% in the field measured. 

F. J. M. 


1036. A. W. VANcE & C. C. SHUMARD,”’ Superregulated power supplies,” 
Electronics, v. 24, no. 12, Dec. 1951, p. 109-113. 


This article describes the power supplies used in the R.C.A. Project 
Typhoon guided missile simulator, built for the U. S. Navy. 
F. J. M. 


1037. G. L. WEIssLER, A. W. Ernarsson & J. D. MCCLELLAND, “‘A semi- 
automatic graphical computer with use in spectroscopy,” Rev. Sct. 
Instruments, v. 23, 1952, p. 209-212. 


The computer is a simple but ingenious device for graphing a prescribed 
function ¢ of a function f(x) where f(x) is given by a graph. The graph of 
y = f(x) appears on a sheet of paper which passes over a line parallel to the 
y axis. Above this line is a transparent plate which moves parallel to the x 
axis. On this plate, one has a graph x = ¢(y) and there is a recorder moving 
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synchronously with the f graph which records the x motion of the plate. The 
operator moves the plate so that the y = f(x) and x = ¢(y) graphs have 
the same ordinate y on the above-mentioned line and hence the desired out- 
put ¢(f(x)) is obtained from the recorder. The article contains details of the 
construction and of the application mentioned in the title. 

F. J. M. 


NOTES 


147.—STABILITY OF DIFFERENCE RELATIONS IN THE SOLUTION OF ORDI- 
NARY DIFFERENTIAL Equations. In a recent communication, J. Topp 
demonstrated the danger of replacing a differential equation, for computa- 
tional purposes, by a difference equation of higher order. H. RUTISHAUSER? 
has since given some general criteria for determining the stability of differ- 
ence approximations to ordinary differential equations. In the present note, 
some standard step-by-step methods of integrating ordinary linear differ- 
ential equations are examined for stability. 

Let a linear differential equation be replaced by a finite difference ap- 
proximation of order  (i.e., one involving » + 1 tabular values). Then the 
nth tabular entry is calculated from 


(1) Yn + + + + = 0, 

where A, Ao, ---, A, are functions of x and of the interval length h. Now 
suppose the errors existing in the entries Ya—pit) Yn-1 ATE En—p, 
€n—p+1, respectively, then the consequent error in y, is where 

(2) én + + Artn_2 + + Aptnp = 0. 

Consider also for convenience that the above errors result entirely from 
errors in the initial values y:, y2, ---, ¥p- Then the general error given by 


equation (2) is 
€n = GAs" + Gedo” + + (n> 


where 4, ds, ---,@, are constants and A, As, ---, A, are the roots of the 
auxiliary equation 
(3) + Ayr?! + Ard? ? + --- +A, = 0. 


The condition for stability is that all the roots of equation (3) lie inside 
or on the unit circle. 
Topp! considered the differential equation 
and its fourth order central difference replacement 
+ (30 12h*)yn_-2 16yn—3 + = 0. 


As h approaches zero, the roots of the corresponding auxiliary equation 


tend to 1,1, 7 — V48, and 7 + V48, the last root quoted being responsible 
for the instability found by Todd. The fourth order backward difference 
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formula, however, is 
(35 + 12h?)yn — 104yn-1 + 1149-2 — S6yn-s + = 0, 


with auxiliary equation roots tending to 1, 1, 11/35, and Vi1/35 as h 
approaches zero. This formula is thus stable at least for sufficiently small h. 
In general, because of the smaller coefficients employed in calculating y, 
from a backward difference formula, the chance of multiplying an error will 
be correspondingly less. Thus in any step-by-step method, stability is more 
likely to result from using backward than central differences. 

The use of backward difference formulae, however, does not ensure sta- 
bility, for consider the equation 


(4) y= 
where y is to be computed for increasing x. Write 
(5) hyo’ = Vyo + V*y0/2 + V¥y0/3 + + V*y0/n, 


leading to the auxiliary equation, 


(6) h+ (1 — 1/d) + (1 — 1/d)?/2 + (1 — 1/d)*/3 + --- 
+ (1 — 1/A)*/n = 0. 


It can be shown that for 4 = 0 the roots of equation (6) other than \ = 1 
have modulus less than unity for n = 6. For nm = 7, there is a pair of con- 
jugate complex roots approximately equal to + 7. For m = 8 there will be 
at least one pair of conjugate roots of modulus greater than unity. (The 
greatest pair is found easily by GRAEFFE’s root-squaring method.) Table I 
demonstrates the instability of the twelfth order backward difference 
formula applied to equation (4). The values of e~* at decimal intervals from 
0 to 1.1 required to start the computations, were taken from five figure 
tables. The theoretical and computed values from 1.2 to 2 are given in rows 
(1) and (2) respectively in table I. 


Table I 


x 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 


(1) .30199 .27253 .24660 .22313 .20190 .18268 .16530 .14957 .13534 
(2) .30148 .27341 .24708 .22244 .20337 .18665 .16017 .14168 .16065 


Next consider the stability of the integration formulae due to ADAMS 
and Mou.ton. Adams’ method is based on the formula 


(7) (91 — Yo)/h = yo’ + 1/2V + + 3/8V*y0' 

+ 251/720V*yo’ + --- 
where a sufficient number of starting values for y, y’ is supposed computed 
by an independent method (e.g., by Taylor series). Consider again the first 
order equation (4). Adams’ formula leads to the auxiliary equation 
(8) FA) =A —1+h{1 + 1/2(1 — 1/d) + 5/12(1 — 1/d)? 

+ 3/8(1 — 1/d)* + 251/720(1 — 1/d)*} = 0, 
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where only fourth differences are retained and it is clear that there is 
stability as 4 tends to zero. Now from equation (8), F(— ©) <0, and 
F(— 1) = — 2+ 551h/45. There is therefore a root of modulus greater 
than unity when h exceeds 90/551, and the method is stable only for suffi- 
ciently small tabular interval. Moreover if higher order differences are re- 
tained, the maximum value of h for which the method is stable is decreased. 

Similar arguments show that Moulton’s method based on the formula 


(9) (yo — y-1)/h = yo! — 1/2V yo" — — 1/24V%y0" 
— 19/720V*yo' — --- 


is also unstable for large values of the tabular interval when differences 
higher than the first are retained. The upper limit on h/ for stability for a 
given number of differences is very much higher than in Adams’ method. 

It has been remarked by RUTISHAUSER? that the error equation, corre- 
sponding to a non-linear differential equation of the form 


= f(x, y™, cee, y*), 


is linear. The above arguments with certain modifications can therefore be 
applied to the stability associated with equations of this form. 
A. R. MItrcHELL 


J. W. Craccs 
Department of Mathematics 
St. Andrews University 
Scotland 
1J. Topp, “Solution of differential equations by recurrence relations,” MTAC, v. 4, 
1950, p. 39-44. 


2H. RuTISHAUSER, “Uber die Instabilitat von Methoden zur gewohnlicher 
Differentialgleichungen,” Zeit. angew. Math. Phys., v. 3, 1952, p. 65-74 


148.—Two NON-ELEMENTARY DEFINITE INTEGRALS. The two integrals 


in question are 
F(x) = f tdt, G(x) = f t‘dt 
0 0 


and are of interest because of the peculiar branching properties of the inte- 
grands and because they lead to series with unusually rapid convergence. 
Integrals of these types have been encountered in some recent studies of 
transients in networks. They can be evaluated numerically as follows. 

As usual, we interpret ¢' as 


et = (t log t)"/n! 


The integral of the general term of this series 


1 
I, = aid, (t log #)"dt 


can be expressed in terms of the complete and incomplete Gamma function 
by means of the transformation 


u = — (n+ 1) logit. 
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In fact 
(- 1) (n + 1)*+17,, =1- T'(n + 1) 
where 
+1) = e~“u"du 
and 


y = — (n + 1) logx. 


However, since available tables! of the incomplete Gamma function are not 
convenient for this application, it was decided to use a direct evaluation of 
the series. 


The integral for I, can be expressed as a sum of m + 1 terms as follows: 
1, = % + — 


Summing over m and collecting the coefficient of each power of log x we 
obtain the double series 


F(x) = In = ¥ (log x)*/r! (— 


r=0 


Similarly 
G(x) = (— = (= 1)"(log x)*/r! 


When x = 1 we have the unusually rapidly converging series 


F(i) 
Gi) 


The following values of F(x) and G(x) were computed to 10D by the above 
double series and checked by numerical integration using tables of fractional 
powers.’ Values have been rounded off to 8D. Most of the calculation was 
made by Mrs. Joan M. Cray. 


x F(x) G(x) 
0 0.0000 0000 0.0000 0000 
A 0.0870 9546 0.1152 6443 
a 0.1625 6731 0.2478 5617 
a 0.2334 0100 0.3890 5018 
A 0.3027 3442 0.5332 8125 
oS 0.3726 1486 0.6763 8768 
6 0.4446 5226 0.8152 2316 
P| 0.5202 8866 0.9474 7031 
8 0.6009 4326 1.0715 0820 
9 0.6881 0902 1.1862 9918 

1.0 0.7834 3051 1.2912 8600 

M. S. CorRINGTON 
Radio Corporation of America 


Camden, N. J. 
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1 PEARSON, Tables of the Incomplete T- Function. London, H. M. Stationery Office. 
(Reissue 1934; London, Biometrika Office, University College. 


) 
2W. GrosNnerR & N. Horreiter, Integraliafel. Erster Teil, Unbestimmte Integrale. 
Vienna, 1949, B: 111, formula 2a. 


3 NBSMTP, Tables of Fractional Powers. New York, Columbia University Press, 1946. 


149.—The Canon DocTRINAE TRIANGVLORVM (1551) of RHETICUS 
(1514-1576). Some facts with reference to this excessively rare publication 
have been given in material about Pitiscus and Rueticus in MTAC, v. 3, 
p. 394, 396, 553-554, 559-560. It is here noted that the only copies known to 
have been preserved were in the Bibliothéque Nationale and British Mu- 
seum. DEMorGAN had a copy in 1845 when he published! a description of 
the work, but this was doubtless in his Library at the University of London, 
destroyed during the recent World War. 

In Catalogue 19, 1952, of the London bookseller E. WEIL, a copy was 
offered for 27 £ 10 s. Mr. WrLL1am D. Moraan, of 1764 St. Anthony Ave., 
St. Paul 4, Minnesota, was so fortunate as to secure this item for adding 
to his already valuable collection (see MTAC, v. 3, p. 562-563). Since Mr. 
Morgan graciously loaned this precious work to me that a microfilm copy 
might be made for the Brown University Library, I take the opportunity 
to add a little to the information already published in MTAC. The com- 
plete title is as follows: Canon Doctrinae Triangvlorum. primom a 
Georgio Ioachimo Rhetico, in lucem editvs, cum privilegio imperiali, Ne quis haec 
intra decennivm, quacunque; forma ac compositione, edere, neve sili vendicere 
aut operibus suis inserere ausit. Lipsiae ex officina Wolphgangi Gunteri. Anno 
M.D. LI. In the title page is an obelisk with a man drawing a diagram on 
the base. 

The back of the title page is blank; then follows a page of Latin verses. 
On the back of this page is the first of 14 pages of 7D tables of the six trigo- 
nometric functions, at interval 10’, arranged for the first time in semi- 
quadrantal display. The degrees are in black, and the minutes and differences 
are in red. This is the first table in which all trigonometric functions are 
brought together. Rheticus was the first to define trigonometric functions 
by means of a right-angled triangle without any reference to a circle. 

Immediately following the tables are 6 pages of dialogue between 
Philomathes, a supposed friend of Rheticus, and Hospes, his pupil. The pupil 
asks what the intention of the book is, and is answered at length. He suggests 
that, perhaps, the intention may be to complete the system of Copernicus, 
by publishing tables from it resembling those then in use. But he is answered 
that Rheticus would rather that Copernicus himself had not done so much 
in this line, as he thereby diminished the geometrical practice of the learner, 
and so on. 

An undated 1580 reprint of the Canon is in the British Museum. 

The copy of the Canon before me has evidently had its pages trimmed; 
but the present size of its pages is 15.8 X 22.5 cm. 

R. C. ARCHIBALD 
Brown University 
Providence, R. I. 


1 DEMorGaN, “On the almost total disappearance of the earliest trigonometrical canon,” 
oe a oe “ 6, 1845, p. 221-228; reprinted with an addition in Phil. Mag., s. 3, v. 26, 
»?P. 
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150.—SQuarRE Root ON THE 602A. The square root setup in the [BM 
602A manual [8th ed., p. 81-84] extracts six digit roots of eight digit base 
numbers at the rate of approximately two roots per minute, using the 
Newton-Raphson scheme with a fixed (six) number of iterations, with 
starting values of 3, 30, 300, or 3000, depending on the size of the base num- 
ber. The 602A must be equipped with division circuits. 

About 75 percent of the 602A’s in use do not have division circuits. 
Square roots can be calculated by many iterative schemes, however, of 
which one of the simplest to program is 


(1) = + — x?) 


where N is the number whose square root is required, and is taken such that 
0 < N <1. The starting value, xo, is taken as N. The iterations are ter- 
minated when 


(2) — xi] <e 


where ¢ is a predetermined small number (say 5 in the 4th decimal place). 
The error is then less than or equal to eN-}. 

The rate of convergence is given by the following inequality where 
N} — x, is the error after the kth iteration: 


|N* — < (Nt — N)(1 + N — 


One notices that convergence is slow for small values of NV. 
For N, a 6 digit number, and e = 10-*, square roots to 4D were obtained 
in the following times. 


N time in seconds 
0.9 25 
0.8 29 
0.7 38 
0.6 50 
0.5 56 
0.4 63 
0.3 75 
0.2 103 
0.1 157 


Starting values, xo, could, of course, be gang punched from a master 
deck, which would make the setup considerably faster. However, the order 
of the deck is then disturbed. If the order of the deck is not critical, the 
remarks below are pertinent; moreover, the entire job (for three digit roots) 
can be done readily from a master deck of 1000 cards, by collating and gang 
punching. 

If division is available and it is possible to order the cards on the numbers 
N, then the formula 


(3) = + N/xi) 
gains speed if the starting value for each card is the result from the previous 
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card. Approximately eight cards per minute can be calculated, with eight 
digit N’s and square roots to 4D. 

The answer from the previous card is taken as xo. The chief problem in 
the wiring is to obtain a starting value for the first card, since the 602A 
performs all the programming wired, in dummy form, before the first card 
is read. If the cards are in descending order on N, the first starting value can 
be taken as N, provided that the dummy programming can be skipped. This 
latter can be accomplished by wiring to ‘“‘read’’ from an early program 
through the normal side of a selector which is then latched for the remainder 
of the run. 

Diagrams of the setups used at the Numerical Analysis Laboratory will 
be published shortly. 


J. HOLLINGSwWorRTH 
F. GRUENBERGER 


Numerical Analysis Laboratory 
University of Wisconsin 
Madison, Wisconsin 


151. A NumericaL Stupy oF A CONJECTURE OF ._KUMMER.—The 
generalization to the cubic case of the well-known (quadratic) Gauss sum 
was first investigated by KUMMER.' He showed that the expression 


(p—1)/2 


(1) Xp=1+2 cos (2xv*/p) 
v=1 


for all p = 1 (mod 3) satisfies the cubic equation 

(2) f(x) = x* — 3px — pA =0 
where A is uniquely determined by the requirements 

(3) 4p = A? + 27B*, A=1 (mod 3). 


Equation (2) clearly has three real roots for each . Kummer classified the 
primes p = 1 (mod 3) according to whether the Kummer sum is the largest, 
middle or smallest root of equation (2). He conjectured that the asymptotic 
frequencies for these classes of p are (in the order above) 3, 4, 3. To check 
this surmise he calculated the first 45 of the x, and found the densities to 
be .5333, .3111, .1556. 

This problem was brought to the attention of the authors by E. ARTIN 
who suggested the desirability of further testing the conjecture since its 
truth wouid have important consequences in algebraic number theory. 

Accordingly the primes p = 1 (mod 3) from 7 through 9,973 were tested. 
We give below a summary of the resulting densities. In this tabulation we 
have arbitrarily divided the primes into six groups of 100 each, designated 
by I, ---, VI and a final group of eleven primes, VII. 
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Number of primes p = 1 (mod 3) such that 


Xp is the Xp is the Xp is the 
largest root middle root smallest root 


Density 


These results would seem to indicate a significant departure from the con- 
jectured densities and a trend toward randomness. 

The method of calculation was this: Each root of (2) lies in one of the 
intervals (— — (— p, + p'), 2p!) as may be seen directly 
from the form of (2) with the help of (3). For each relevant p the expression 
(1) for x, was evaluated, its sign was determined and its square compared 
to p. This determined in which of the three intervals just described the x, 
lies. To check that x, was indeed a solution of (2), (3) and to determine the 
precision of the evaluation the expression 


(4) (Xp) 


was then calculated. This latter check was performed by first transforming 
the x, into decimal form for tabulation and then retransforming these results 
back into binary form before evaluating the expression (4). In this manner 
both the calculation proper and the conversion to decimal form of the 
results were checked. 

The trigonometric expressions appearing in (1) were evaluated by power 
series. Each angle was reduced mod 27 and then mod z until it lay between 
— 2/2 and + 2/2. Then the cosine of } of this angle was calculated keeping 
five terms in the series expansion. The ‘‘double-angle’”’ formula for cosines 
was then used twice to obtain the desired cosine. 

The calculation involved about 15 million multiplications counting the 
checking mentioned above. The values of p were introduced in blocks of 200. 
The entire calculation was carried out twice to ensure reliability. The 
authors are indebted to Mrs. ATLE SELBERG who programmed and coded 
the calculation. 

J. von NEUMANN 


Institute for Advanced Study H. H. GOLDsTINE 


Princeton, New Jersey 


1E, E. Kummer, “De residuis cubicis disquisitiones nonnullae analyticae,” Jn. f. d. 
reine u. angew. Math., v. 32, 1846, p. 341-365. 


CORRIGENDA 


. 6, p. 262, insert Emch, G. F. 247. 

. 6, p. 265, under Myers insert 54. 

. 6, p. 268, under Yowell insert 254. 

. 7, p. 31, 1. 11 of MTE 218; for — log p read log p. 
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I 54 28 18 : | 
I 41 38 21 
Il 46 33 21 
: IV 39 32 29 
V 43 29 28 
VII 5 3 3 
Total 272 201 138 
| A452 3290 2258 
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CLASSIFICATION OF TABLES 


Arithmetical Tables. Mathematical Constants 
Powers 

Logarithms 

Circular Functions 

Hyperbolic and Exponential Functions 
Theory of Numbers 

Higher Algebra 

Numerical Solution of Equations 
Finite Differences, Interpolation 
Summation of Series 

Statistics 

Higher Mathematical Functions 

. Integrals 


Interest and Investment 


Physics, Geophysics, Crystallography 

Chemistry 

Navigation 

Aerodynamics, Hydrodynamics, Ballistics 
Calculating Machines and Mechanical Computation 


A. 
B 
Ags ol 
D. 
E. 
F, 
G. 
K, 
a 
M 
P. Engineering 
Q. Astronomy 
R. Geodesy 
S. 
U. 
Vv. 
Z. 
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